L'&UUUUU 


00 

00 

<M 

(M 

hs 


< 

I 

D 

< 


David  W.  Taylor  Naval  Ship  Research  and  Development  Center 

Bethesda,  MD  20084-5000 


DTNSRDC- 86/054  September  1986 

Ship  Performance  Department 
Research  and  Development  Report 


Far-Field  Features 
of  the  Kelvin  Wake 

by 

A.  Barnell 
F.  Noblesse 


THIS  DOCUMENT  IS  BEST 

QUALITY  AVAILABLE.  THE  COPY 

* 

FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 


MAJOR  DTNSRDC  TECHNICAL  COMPONENTS 


COMMANDER  00 

- ! - 

TECHNICAL  DIRECTOR  01 


OFFICER  IN  CHARGE 
CARDEROCK 

05 

SHIP  SYSTEMS  INTEGRATION 
DEPARTMENT 

12 

SHIP  PERFORMANCE 
DEPARTMENT 

15 

AVIATION  AND  SURFACE 
EFFECTS  DEPARTMENT 

16 

STRUCTURES  DEPARTMENT 

17 

COMPUTATION,  MATHEMATICS 
&  LOGISTICS  DEPARTMENT 

18 

SHIP  ACOUSTICS  DEPARTMENT 

19 

CENTRAL  INSTRUMENTATION 
DEPARTMENT 

29 

OFFICER  IN  CHARGE 
04  ANNAPOLIS 


t  PROPULSION  AND  AUXILIARY 
27  SYSTEMS  DEPARTMENT 


SHIP  MATERIALS  ENGINEERING 
28  DEPARTMENT 


DESTRUCTION  NOTICE  —  For  classified  documents,  follow  the  procedures  in  DOD 
5220.22M,  Industrial  Security  Manual,  Section  11-9,  or  DOD  5200.1-R,  Information 
Security  Program  Regulation,  Chapter  IX.  For  unclassified,  limited  documents, 
destroy  by  any  method  that  will  prevent  disclosure  of  contents  or  reconstruction  of 
the  document. 


-.X. 


.  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE 


la  REPORT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


Alia  vr 


REPORT  DOCUMENTATION  PACE 


1b.  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION/ DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
DTNSRDC-86/ 054 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION  1 6b  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 


David  Taylor  Naval  Ship 
R&D  Center 


6c  ADDRESS  (City.  State,  and  ZIP  Code) 
Bethesda,  MD  20084-5000 


(If  applicable) 
Code  1542 


7b.  ADDRESS  (City.  State,  and  ZIP  Code) 


8a  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


fit  ADDRESS  (City,  Hate,  and  ZIP  Code) 


II  HUE  (Include  Security  Classification) 

V  Far-Field  Features  of  the  Kelvin  Wake 


17  personal  AUTHOR(S) 

Barnell,  A.  and  Noblesse,  F. 


to  SOURCE  OF  FUNDING  NUMBERS 


Ha  TYPE  OF  REPORT 

Final 


13b  TIME  COVERED 
FROM  TO 


PROJECT 

TASK 

NO 

NO, 

RS43421 

RS43421 

WORK  UNIT 
ACCESSION  NO 

DN178067 


14  DATE  OF  REPORT  (Year.  Month.  Day)  15  PAGE  COUNT 

1986,  September  33 


COSATI  CODES _ I  18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

roup  I  sub-group  I  Ship  Waves  Divergent  Waves 

Kelvin  Wake 


'9  ABSIR-W.T  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  classical  Kelvin  wake  trailing  a  ship  advancing  at  constant  speed  in  calm  water 
is  studied.  In  particular,  asymptotic  expressions  for  determining  the  far/field  wake,  at 
large  distances  behind  the  ship,  are  investigated  numerically.  This  analysis^-dn  which8-' 
rfaee-^rens'ion-and-nenWnea r it ies ■  indicates  that  the  divergent  waves  o  a 

surface  ship  generally  are  infinitely  steep  at  the  track  of  the  ship,  even  though  their 
amplitude  vanishes  there;  this  result  is  theoretically  possible  because  the  wavelength  of 
the  divergent  waves  vanishes  at  the  track  of  the  ship^  Inasmuch  as  infinitelyfsteep  water 
waves  cannot  exist  in  reality,  the  linear  noisurface/tension  analysis  presented  in  this 
study  suggests  that  no  divergent  waves  can  exist  within  a  certain  region  in  the  vicinity 
of  the  track  of  the  ship,  and  that  the  Kelvin  wake  contains  threie  distinct  regions. 

,^(4^  an  inner  region  where  only  transverse  waves  can  exist}  (wf  an  outer  region  where  both 
transverse  and  divergent  waves  are  present/  and  (• irxfan  intermediate  region  at  the  bound¬ 
ary  between  the  inner  and  outer  regions  where  short  steep  divergent  waves,  as  well  as  -  i 

(Continued) 


20  DiSIRIBUIlON  /  AVAILABILITY  OF  ABSTRACT 

09 UNCI ASSlFIEO'UNLlMITEO  □  SAME  AS  RPT  DOTIC  USERS 


77 a  NAME  OF  RESPONSIBLE  INDIVIDUAL  ZZb  TELtrH 

Alexander  Barnell  ^(301) 

DD  FORM  1 473, 84  mar  83 APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED.  _ 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOl 

(301)  227-1865  Code  1542 


SECURITY  CLASSlf  ICA I  ION  Ol  THIS  PAGE_ 
UNCLASSIFIED 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS'PAGE 


1 


(Block  19  continued) 

transverse  waves,  can  be  found.  Numerical  results  for  a  simple  bow  form  show  that  the 
inner  region  is  quite  narrow,  and  that  the  wavelength  of  the  divergent  waves  at  the 
boundary  of  the  inner  region  is  of  the  order  of  1  to  20  cm,  depending  on  the  speed  of 
the  ship.  These  results  appear  to  be  consistent  with  the  narrow  V-wakes  observed  in 
some  SAR  images  of  ship  wakes . The  analysis  presented  in  this  study  may  thus  provide  a 
partial  explanation  for  these  observations.  However, 'the  small  wavelengths  mentioned 
above  indicate  that  the  divergent  waves  in  the  vicinity  of  the  track  of  the  ship  are 
likely  to  be  significantly  altered  by  surface  tension,  which  -should  "thetef or^  be  included 


in  a  more  realistic  analysis. 


(A4,U..S 


Accession  For 


NTI6  GRAfcl 
DIIC  TAB 
Unannounced 
Justification. 


¥ 

□ 


Distribution/ 


Availability  Codes 
Avail  and/or 
Special 


4 


UNCLASSIFIED  _ _ 


lU  CljlM 1  v«  A* 


FIGURES  (Continued) 


Page 

8.  Waterlines  and  framelines  of  the  simple  ship  bow  form 

considered  for  numerical  applications... . .  8 

9a.  Real  and  imaginary  parts  of  the  functions  K{{(t)  and  K()(t)  for  a 
simple  ship  bow  form  with  8  =  12°  (top)  and  3°  (bottom), 

Y  =  0  and  F  =  0.5 . . . . .  11 

9b.  Real  and  imaginary  parts  of  the  functions  Kj^(t)  and  Ko(t)  for  a 

simple  ship  bow  form  with  0=1°  (top)  and  20*  (bottom), 

Y  3  0  and  F  3  0.5 . . . . . . .  11 

10.  Real  and  imaginary  parts  of  the  functions  K^tt)  and  ^0^0  for  a 

simple  ship  bow  form  with  Y  3  0  and  45°,  8  *  12°  and  F  3  0.5 .  12 

11.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  %(t)  and 

Ko(t)  for  a  simple  ship  bow  form  with  0  *  12°,  Y  *  45°  and  F  3  0.5 .  12 

12a.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  K±a£/ 

(+0+1)1^2  for  k  =  0  (left)  and  1  (right),  K  3  KM  (Michell 
thin-ship  approximation)  and  a  simple  ship  bow  form  with 

0  3  12°,  Y  3  45°  and  F  3  0.5 .  13 

+  . 

12b.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  K+a£/ 

(+Q+')1/2  for  k  =  2  (left)  and  3  (right),  K  3  KM  (Michell 
thin-ship  approximation)  and  a  simple  ship  bow  form  with 

0  =12°,  Y  3  45°  and  F  3  0.5 .  H 

4. 

13a.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  K+aj^/ 

(+0V)1/2  for  k  3  0  (left)  and  1  (right),  K  3  Kq  (zeroth-order 

slender-ship  approximation)  and  a  simple  ship  bow  form  with 

0  =  12°,  Y  3  45°  and  F  3  0.5 .  13 


CONTENTS 


Pa 

ABSTRACT . . . . . . . . 

INTRODUCTION . , . . . 

APPROACH . . . . . . 

ASYMPTOTIC  EVALUATION  OF  THE  KELVIN  WAKE . . . . . . 

A  SIMPLE  TEST  CASE:  THE  FAR-FIELD  WAVE-AMPLITUDE  FUNCTION . .  .  . . 


A  SIMPLE  TEST  CASE:  THE  FAR-FIELD  KELVIN  WAKE... . . .  1 

SUMMARY  OF  RESULTS  AND  CONCLUSIONS... . 1 

ACKNOWLEDGMENTS . 1 

REFERENCES..... .  1 

FIGURES 


1.  Real  parts  of  the  functions  E+(t;x,ri)  and  E_(t;x,ot)  for  x  *  -50  and 

a  ■  .1,  .2,  1/23/2,  ,4,  and  .5.... . 

2.  Real  and  imaginary  parts  of  the  functions  (l+t2)K(t)  and  t(l+t2)K(t) 

for  a  simple  ship  bow  shape  in  the  zeroth-order  slender-ship 
approximation . . . . 

3a.  The  phase  function  0(t;a)  for  0  <  t  <  6  and  several  values  of  a  .... 

3b.  The  function  9'(t;a)  for  0  <  t  <  6  and  several  values  of  ot  . 

3c.  The  function  9"(t;a)  for  0  <  t  <  3  and  several  values  of  a  . 

4.  The  functions  tg(cx),  t+(o)  and  t_(ci)...... . . . 

5.  The  functions  tg(a),  t±(a),  Og(a),  Oq*  '  ^ <x) »  ®±*  (a) 

and  0+''(ct)  for  0  <  a  <  0.5 . . . . . 

6.  The  classical  Kelvin  ship  wave  pattern . . . 

7.  The  wavelengths  A+(a )  and  propagation  angles  fe(oi)  of  the  transverse 

and  divergent  waves  in  the  Kelvin  wake . . 


iii 


FIGURES  (Continued) 


13b.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  K±a^/ 
(TQ+')*/2  for  k  ■  2  (left)  and  3  (right),  K  *  Kq  (zeroth-order 
slender-ship  approximation)  and  a  simple  ship  bow  form  with 


g  -  12°,  Y  -  45°  and  F  -  0.5. 


14.  Real  (top)  and  imaginary  (bottom)  parts  of  the  functions  K±a4/ 
(^©^ » )l/2  £or  (left)  and  Kq  (right)  and  a  simple  ship 


bow  form  with  g  ■  12°,  Y  *  45°  and  F  -  0.5. 


15.  Amplitude  of  the  transverse  (-)  and  divergent  (+)  waves  in  the  Kelvin 
wake  for  a  simple  ship  bow  form,  with  g  “  12°,  Y  m  45°  and  F  *  0.5, 
in  the  Michell  thin-ship  approximation  (left)  and  the  zeroth-order 


slender-ship  approximation  (right), 


16.  Steepness  of  the  transverse  (~)  and  divergent  (+)  waves  in  the  Kelvin 
wake  for  a  simple  ship  bow  form,  with  g  *  12°,  Y  »  45°  and  F  -  0.5, 
in  the  Michell  thin-ship  approximation  (left)  and  the  zeroth-order 


slender-ship  approximation  (right). 


17.  Upper  bound  for  the  steepness  of  the  divergent  wave3  in  the  Kelvin 

wake  for  a  simple  ship  bow  form,  with  g  *  12°,  Y  ■  45°  in  the  Michell 
thin-ship  approximation  (left)  and  the  zeroth-order  slender-ship 


approximation  (right), 


18.  Steepness  of  the  divergent  waves  in  the  Kelvin  wake  for  a  simple 
ship  bow  form,  with  g  =  12° ,  Y  =  0  (top)  and  45°  (bottom),  and 
F  =  0.3  (left)  and  0.8  (right)  in  the  zeroth-order 


slender-ship  approximation. 


v 


FIGURES  (Continued) 


Page 

19.  Amplitude  of  the  peak  value  of  the  steepness  of  the  divergent  waves 
in  the  Kelvin  wake  for  a  simple  ship  bow  form,  with  (3  *  12°  and 

Y  3  45°  in  the  zeroth-order  slender-ship  approximation....... . .  16 

20.  The  Kelvin  cusp  line,  the  line  along  which  the  amplitude  of  the 
divergent  waves  exhibits  a  peak,  and  the-  three  lines  along 
which  the  steepness  of  the  divergent  waves  is  equal  to 

1/20,  1/15  and  1/7  (chain  line  close  to  the  track  of 
the  ship)  for  a  simple  ship  bow  form,  with  0  *  12® 
and  Y  3  45®,  in  the  zeroth-order  slender-ship 

approximation. . . . . . . . . . . . .  16 

21.  Lines  along  which  the  steepness  of  the  divergent  waves  is  equal  to 
1/20,  1/15  and  1/7  for  a  simple  ship  bow  form,  with  0  *  12®  and 

Y  3  45®,  in  the  Michell  thin-ship  approximation  ( - )  and  the 

zeroth-order  slender-ship  approximation  ( - r) . . .  17 

22.  Lines  along  which  the  steepness  of  the  divergent  waves  is  equal  to 


1/15  for  a  simple  ship  bow  form  with  0*  8°,  12®,  16®  and 
Y  3  45°  (top)  and  Y  =  0,  25°,  45°  and  0  =  12°  (bottom) 
in  the  zeroth-order  slender-ship  approximation... . 


17 


FAR-FIELD  FEATURES  OF  THE  KELVIN  WAKE 


A.  BotmN  and  F.  NoMrae 

David  W.  Taylor  Ni«il  Skip  Researtk  and 'Development  Center 
Betfceada,  Maryland  20M4 


ABSTRACT 

The  classical  Kelvin  wake  trailing  a  ship  advancing 
at  constant  speed  in  calm  water  is  studied.  In  particular, 
asymptotic  expressions  for  determining  the  far-field  wake, 
at  large  distances  behind  the  ship,  are  investigated 
numerically.  This  analysis,  in  which  surface  tension  and 
nonlinearities  are  ignored,  indicates  that  the  divergent 
waves  of  a  surface  ship  generally  are  infinitely  steep  at  the 
track  of  the  ship, .even  though  their  amplitude  vanishes 
there;  this  result  is  t heoret ically 1  possible  because  the 
wavelength  of  the  divergent  waves  vanishes  at  the  track  of 
the  ship.  Inasmuch  as  infinitely-steep  water  waves  cannot 
exist  in  reality,  the  linear  no-surface-tension  analysis 
presented  in  this  study  suggests  that  no  divergent  waves 
can  exist  within  a  certain  region  in  the  vicinity  of  the 
track  of  the  ship,  and  that  the  Kelvin  wake  contains  three 
distinct  regions:  (i)  an  inner  region  where  only  transverse 
waves  can  exist,  (ii)  an  outer  region  where  both  transverse 
and  divergent  waves  are  present,  and  (iii)  an  intermediate 
region  at  the  boundary  between  the  inner  and  outer 
regions  where  short  steep  divergent  waves,  as  well  as 
transverse  waves,  can  be  found.  Numerical  results  for  a 
simple  bow  form  show  that  the  inner  region  is  quite 
narrow,  and  that  the  wavelength  of  the  divergent  waves  at 
the  boundary  of  the  inner  region  is  of  the  order  of  I  to 
20  cm,  depending  on  the  speed  of  the  ship.  These  results 
appear  to  be  consistent  with  the  narrow  V-wakes  observed 
in  some  SAR  images  of  ship  wakes.  The  analysis 
presented  in  this  study  may  thus  provide  a  partial 
explanation  for  these  observations.  However,  the  small’ 
wavelengths  mentioned  above  indicate  that  the  divergent 
waves  in  the  vicinity  of  the  track  of  the  ship  are  likely  to 
be  significantly  altered  by  surface  tension,  which  should 
therefore  be  included  in  a  more  realistic  analysis. 


INTRODUCTION 

It  has  been  observed,  e.g.  [1,2,3],  that  images  of 
ship  wakes  taken  by  the  SEASAT  SAR  (Synthetic 
Aperture  Radar)  and  other  airborne  SAR  systems 
sometimes  reveal  long,  narrow  wakes  extending  back 
several  miles  behind  ships.  According  to  Swanson  [3], 
“there  appears  a  narrow,  dark  band  along  the  centerline 
of  the  ship  wake,  widening  slightly  with  distance  away 


from  the  ship.  Often  accompanying  these  dark  regions,  at 
a  slightly  larger  angle,  are  bright  lines,  which  show  up 
especially  well  on  X-band  SAR  images".  The  apparent 
included  angle  of  this  narrow  wake  is  much  smaller  that 
that  corresponding  to  the  Kelvin  cusp  line. 

Several  theoretical  explanations  of  the  features 
observed  in  SAR  images  of  ship  wakes  have  been 
proposed  [3, 4, 5, 6].  The  proposed  explanations  include 
interactions  between  the  cross  currents  created  by  a  ship  in 
its  wake  and  surface  gravity  waves  |I,6)  and  the  occurence 
of  a  sharp  peak  in  the  amplitude  of  the  divergent  waves 
in  the  Kelvin  wake  for  a  ship  form  having  a  large  flare 
angle  [5],  However,  there  currently  appears  to  be  no 
convincing  electrodynamic  ot  hydrodynamic  explanation  of 
the  observations.  It  thus  seems  necessary  to  investigate  in 
detail  each  one  cf  the  proposed  possible  explanations  in 
order  to  determine  whether  one,  or  possibly  a  combination 
of  several,  of  them  does  in  fact  provide  a  satisfactory 
theoretical  interpretation  of  the  SAR  images  of  ship 
wakes. 

The  classical  Kelvin  wake  trailing  a  ship  advancing 
at  constant  speed  in  calm  water  is  examined  in  this  study. 
Viscosity,  surface  tension  and  nonlinearities  are  ignored  in 
this  preliminary  analysis,  and  irrotation&l  flow  is  assumed. 
The  analysis  and  related  numerical  results  for  a  simple 
ship  bow  form  have  led  to  several  conclusions  and 
recommendations  for  further  studies  which  are  presented 
at  the  end  of  the  paper.  Only  two  conclusions  that  may 
be  related  to  some  of  the  features  observed  in  SAR 
images  of  ship  wakes  are  discussed  here. 

it  was  previously  found  by  Scragg  [5]  that,  for  a 
ship  bow  form  with  a  large  flare  angle,  the  zeroth-order 
slender-ship  approximation  to  the  far-field  wave-amplitude 
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function  liven  in  Noblesse  [7]  predicts  a  sharp  peak  in  the 
value  of  the  amplitude  of  the  divergent  wives  at  an  angle 
from  the  track  of  the  ship  equal  to  approximately  half  the 
bow  entrance  angle.  This  finding  of  Scragg  has  been 
verified  in  this  study,  as  may  be  seen  from  figures  IS,  16 
and  18.  Furthermore,  the  magnitude  of  the  wave- 
amplitude  peak  has  been  found  to  increase  very  rapidly  as 
the  Froude  number  decreases  below  a  certain  threshold 
value.  This  peak  thus  appears  to  be  a  large-flare  low- 
Froude-number  feature.  For  the  simple  ship  bow  form 
considered  in  this  study,  the  peak  in  the  amplitude  of  the 
divergent  waves  in  the  Kelvin  wake  occurs  along  a  line  at 
6°  from  the  track  of  the  ship.  It  may  be  found  from 
equation  (30a)  that  the  wavelength  of  the  divergent  waves 
along  this  line  varies  between  0.7  m  and  4.3  m  for  a  ship 
with  speed  varying  between  10  knots  and  2S  knots, 
respectively. 

The  other  conclusion  of  this  study  that  may  be 
related  to  the  narow  V-wakes  observed  in  some  SAR 
images  of  ship  wakes  is  the  result  that  the  divergent  waves 
of  a  surface  .hip  are  infinitely  steep  at  the  track  of  the 
ship,  even  though  their  amplitude  vanishes  there.  This 
result  is  theoretically  possible  because  the  wavelength  of 
the  divergent  waves  vanishes  at  the  track  of  the  ship.  A 
similar  result  was  previously  obtained  by  Sharma  18)  who 
found  that  the  Michell  thin-ship  approximation  for  a  thin 
and  deep  strut-like  ship  form  predicted  infinite  slopes  for 
the  divergent  waves  at  the  track  of  the  ship.  Inasmuch  as 
infinitely-steep  water  waves  cannot  exist  in  reality,  the 
foregoing  result  suggests  that  no  divergent  waves  can  exist 
within  a  certain  region  in  the  vicinity  of  the  track  of  the 
sh'p,  and  that  the  Kelvin  wake  contains  three  distinct 
regions:  (i)  an  inner  region  where  only  transverse  waves 
can  exist,  (ii)  an  outer  region  where  both  transverse  and 
divergent  waves  are  present,  and  (iii)  an  intermediate 
region  at  the  boundary  between  the  inner  and  outer 
regions  where  short  steep  divergent  waves,  as  well  as 
transverse  waves,  can  be  found. 

Numerical  results  for  a  simple  bow  form  show  that 
the  “no-divergem-wave"  inner  region  is  quite  narrow,  as 
may  be  seen  from  figure  20  showing  the  Kelvin  cusp  line 
(angle  y  19°28’).  the  line  along  which  the  amplitude  of 
the  divergent  waves  exhibits  a  peak  (angle  y  6°),  and  the 
three  lines  along  which  the  steepness  of  the  divergent 
waves  is  equal  to  1/20,  1/15  and  1/7  (chain  line  close  to 
the  track  of  the  ship).  The  latter  three  lines,  along  which 
the  divergent  waves  are  steep,  lie  much  closer  to  the  track 


of  the  ship  than  the  line  along  which  the  steepness  of  the 
divergent  waves  exhibits  a  peak.  Figures  21  and  22  show 
considerable  variations  among  the  several  “steep-divergent- 
waves"  lines  that  are  represented  in  these  figures. 
Nevertheless,  these  lines  may  be  seen  to  correspond  to 
values  of  Y/(-X)  equal  to  about  10-2  to  2.10-2.  For 
ship  speeds  varying  between  10  knots  and  25  knots,  the 
wavelength  of  the  divergent  waves  corresponding  to  the 
foregoing  values  of  Y/(-X)  may  be  shown  to  vary 
between  0.7  cm  and  4.2  cm  for  Y/(-X)  =  10~2,  and 
2.7  cm  and  17  cm  for  Y/(-X)  =  2.10'2.  These 
wavelengths  are  consistent  with  the  wavelengths  of  the 
radar  pulses  used  in  SAR  imaging,  so  that  the  foregoing 
results  may  provide  a  partial  hydrodynamic  explanation 
for  the  narrow  V-wakes  observed  in  these  images. 

The  tentative  nature  of  this  explanation  must 
however  be  stressed.  Indeed,  the  foregoing  results  are 
based  on  an  analysis  in  which  surface  tension  and 
nonlinearities  have  been  neglected.  Inasmuch  as  this  linear 
no-surface-tension  analysis  predicts  extremely  short  and 
steep  waves  in  the.vicinity  of  the  track  of  the  ship,  it  is 
evident  that  both  surface  tension  and  nonlinear  effects  are 
liable  to  be  significant.  In  particular,  the  short 
wavelengths  found  along  the  steep-divergent-wave  lines 
determined  in  this  study,  and  the  brief  description  of  the 
effects  of  surface  tension  upon  the  Kelvin  wake  given  in 
Sharma  [8),  Lamb  (9,  pp.  468-470)  and  Wehausen  and 
Laitone  (10,  pp.  636-637)  indicate  that  the  system  of 
divergent  waves  in  the  vicinity  of  the  track  of  the  ship  js 
likely  to  be  profoundly  affected  by  surface  tension. 

Effects  of  surface  tension  upon  the  Kelvin  wake  will  be 
investigated  in  a  sequel  to  the  present  study. 

APPROACH 

This  study  considers  the  steady  potential  flow  due 
to  a  ship  advancing  with  constant  speed  in  calm  water  of 
infinite  depth  and  lateral  extent.  The  far-field  Kelvin 
wake,  which  is  of  primary  interest  here,  may  be 
conveniently  analyzed  in  terms  of  the  nondimensional  far- 
field  coordinates  T  =  Xg/U2,  velocity  potential  |  =  ♦g/ll3 
and  velocity  vector  Vxj  =  Vx<h/U,  where  g  is  the 
gravitational  acceleration  and  U  is  the  speed  of  advance  of 
the  ship,  X  and  <t>  represent  the  dimensional  coordinates 
and  velocity  potential,  respectively,  and  Vx  and  Vx  are  the 
nondimensional  and  dimensional  differential  operators 
Vx  =  (d/dx,d/dy,d/dz)  and  ?x  =  (d/dX,d/d\,d/dZ). 
The  mean  free  surface  is  taken  as  the  plane  z  =  0,  with 
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the  z  axis  pointing  upwards,  and  the  x  axis  is  chosen  in 
the  ship  centerplane  and  pointing  towards  the  bow.  The 
origin  of  the  system  of  coordinates  is  placed  within  the 
ship.  The  Froude  number  is  denoted  by  F  =  U/(gL)l/2, 
where  L  is  the  length  of  the  ship. 

Equation  (32)  in  (7)  yields  the  following  expression 
for  the  velocity  potential  associated  with  the  Kelvin  wake 
behind  the  ship 


n+(x)  =  !m  /  [E+(tix)  +  E_(t;xj]K{t)  dt, 


0) 


where  E^tix)  is  the  exponential  function 
E±(tix)  =  exp(z(l+t2)  +  i(x±ytXl+t2)l/2!,  (2) 

and  K(i)  is  the  far-field  wave-amplitude  function,  which 
depends  on'the  hull  shape  and  the  Froude  number. 
Assuming  that  differentiation  under  the  integral  sign  is 
permitted  in  equation  (1),  we  may  obtain 


=  Im 


.  /°o  E.+E.  (l+t2)1/2 

'/>  E  -E_  K(°  t(l+t2)''2  ** 


(3a) 
(3b) 

The  nondimensional  elevation  e  =  Eg/U2  of  the 
free  surface  at  a  sufficiently-large  distance  behind  the  ship, 
such  that  nonlinearities  may  be  neglected,  is  given  by 
e(x.y)  =  a+(x,y,0)/3x.  (4) 

The  slopes  of  the  free  surface  in  the  directions  parallel 
and  perpendicular  to  the  ship  course  then  are 
3e(x,y)/0x  =  324(x,y,0)/3x2  (4a) 

3e(x,y)/3y  =  32+(x,y,0)/3x3y.  (4b) 

If  differentiation  under  the  integral  sign  in  equation. (3a)  is 
permitted,  we  have 


=  -Im/* 
•'0 


|E++E_ 

!e+-e. 


K(t) 


1+t2 

t(l+t2) 


dt. 


(5a) 

(5b) 


The  vertical  velocity  is  given  by 

4,  =  (6) 

In  this  study,  we  are  mostly  interested  in  the  value 
of  the  several  flow  variables  defined  above  at  the  mean 
free  surface  and  behind  the  ship,  so  that  we  have  z  -  0 
and  x  <  0.  Expression  (2)  for  the  exponential  function  E  ± 
then  becomes 

E±  =  explix(lTat)(l+t2)l/2l,  (7) 

where  a  is  defined  as 


a  =  -y/x.  (8) 

For  a  ship  with  port-  and  starboard-symmetry,  as  is 
considered  here,  the  Kelvin  wake  is  symmetric  about  the 
ship  track  y  =  0.  We  may  then  restrict  the  analysis  of  the 
Kelvin  wake  to  the  domain  y  >  0  and  x  <  0,  and  assume 
a>  0. 


Two  difficulties  associated  with  the  foregoing 
approach  for  numerically  determining  the  potential,  its 


gradient,  and  the  free-surface  slopes  in  the  Kelvin  wake 
are  readily  apparent  and  should  be  noted  here.  A  first 
numerical  difficulty  stems  from  the  oscillations  of  the 
exponential  function. E+  given  by  equation  (7),  which  arc 
very  rapid  for  large  values  of  |x|.  We  have  x  =  Xg/U2 
=  (X/L)/E2.  For  a  typical  value  of  the  Froude  number 
equal  to  0.2,  say,  we  thus  have  x  =  25  X/L;  we  then 
have  x  =  -  250  at  10  ship  lengths  behind  the  ship,  and 
much  larger  values  of  |x|  must  be  considered  at  greater 
distances  behind  the  ship  or/and  for  smaller  values  of  the 
Froude  number.  Even  for  the  comparatively-moderate 
value  of  x  equal  to  -  50,  figure  1  shows  that  the 
functions  E+  and  E_  oscillate  quite  rapidly.  More 
precisely,  figure  1  depicts  the  real  parts  of  the  functions 
E+(t;x,«)  and  E_(t;x,o)  for  x  =  -50;  a  =  .1,  .2, 
1/2J/2,  .4  and  .5;  and  for  0  <  t  <  7  and  0  <  t  <  3  on 
the  left  and  right  sides,  respectively.  Figure  1  also 
indicates  that  the  behavior  of  the  function  E+(t;x,a) 
strongly  depends  on  the-value  of  a. 


*+  (««-»)  e_ 


E_(t;x,a)  for  x  =  -50  and  a  =  .1,  ,2,  l/23/2,  .4  and  .5 


A  second,  more  basic,  difficulty  is  associated  with 
the  differentiation  under  the  integral  sign  which  was  used 
for  obtaining  expressions  (3a, b)  and  (5a, b)  from 
expression  (1).  For  a  fully-submerged  body,  the  far-field 
wave-amplitude  function  K(t)  is  exponentially-smal!  as  t  — 
oo,  so  that  the  operation  of  differentiating  under  the 
integral  sign  in  expression  (1)  can  be  continued  ■•'•’-“finitely 
in  principle.  Differentiation  under  the  integral  sign  ikewise 
is  justified  if  z  <  0.  However,  the  operation  must  be 
justified  in  the  limiting  case  z  =  0  for  a  surface  ship. 
Clearly,  the  operation  may  not  be  justified  in  principle,  or 
feasible  in  practice,  if  the  far-field  wave-amplitude 
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function  K(t)  does  not  vanish  sufficiently  rapidly  as  t  -» 
<*>.  Precise  information  about  the  asymptotic  behavior  of 
the  function  K(t)  as  t  -►  «  is  required  in  this  respect. 
Figure  2  depicts  the  real  and  imaginary  parts  of  the 
functions  (1  +t2)K(t)  and  t(l  +  t2)K(t)  appearing  in  the 
integrands  of  the  integrals  (5a)  and  (5b),  respectively,  for 
0  <  t  <  16  and  for  a  simple  ship-bow  shape  which  is 
considered  further  on  in  this  study.  Differentiation  under 
the  integral  sign  in  expression  (3a)  is  clearly  not  justified 
in  the  case  corresponding  to  figure  2. 
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Fig.  2  —  Real  and  Imaginary  Parts  of  the  Functions 
(1  +t2)K(t)  and  t(l  +t2)K(t)  for  a  Simple  Ship  Bow  Shape 
in  the  Zeroth-Order  Slender-Ship  Approximation 


ASYMPTOTIC  EVALUATION  OF  THE 
KELVIN  WAKE 

For  large  values  of  |x|  and  z  =  0,  analytical 
approximations  to  the  integrals  (1),  (3a, b)  and  (5a, b)  can 
be  obtained  by  taking  advantage  of  the  rapid  oscillations 
of  the  exponential  functions  E±  defined  by  equation  (7). 
These  functions  may  be  expressed  in  the  form 
E  =  exp[ix0(t;a)],  (9) 

where  0(t;a)  is  defined  as 

0(t;a)  =  (1-atXl+t2)1'2.  (10) 

The  functions  E+(t;x,a)  and  E_(t;x,a),  where  a  >  0, 
correspond  to  the  function  E(t;x,a)  with  <*  >  0  and  a  <  0, 
respectively.  In  the  limiting  case  a  =  0,  we  have  E+  = 
exp[ix(l  +  t2)l/2]  =  E_.  The  derivatives  of  the  phase- 
function  0(t;x,a)  with  respect  to  t  are  given  by 
0'(t;a)  =  — (a  — t  +  2at2)/(l+t2)1/2,  (11a) 

0"(t;a)  =  ( 1  -  3at  -  2ot3)/(  1  + 12)3/2,  (11b) 

0"'(t;a)  =  -3(a  +  t)/(l  +  t2)5/2.  (11c) 

The  phase  0(t;a)  and  its  first  and  second  derivatives 


0'(t;«)  and  0"(t;a)  are  depicted  in  figures  3a,  b,  c, 
respectively. 


Fig.  3a  —  The  Phase  Function  0(t;a)  for  0  <  t  <  6  and 
Several  Values  of  a 


Fig.  3b  —  The  Function  8'( t;a)  for  0  <  t  <  6  and 
Several  Values  of  a 


Fig.  ic  —  The  Function  0"(t;a)  for  0  <  t  <  3  and 
Several  Values  of  a 
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It  may  be  seen  from  figure  3b  that  the  first 
derivative  9'  does  not  vanish  except  if  0  <  a  <  l/23/2,  • 
for  which  9'  vanishes  for  2  values  of  t.  Equation  (11a) 
shows  that  the  2  points  white  9'  --0,  that  is  where  the 
phase  9  is  stationary,  are  given  by 
t±(«)  =  [l±(l-8a2)1/2)/4*.  (12) 

We  have  0  <  t_  <  t+  <  »,  with  t_  =  0  and  tH  =  °°  for 
a  =  0,  and  t_  =  1/2I/2  =  t+  for  a  =  l/23/2  The  2 
points  of  stationary  phase  t_  and  t+  are  apparent  on  the 
left  side  of  figure  1  for  a  =  .1,  .2  and  1/2372.  Figure  3c 
shows  that  the  second  derivative  9"  vanishes  for  one 
value  of  t,  say  tg,  if  a  >  0.  Equation  (lib)  yields 
t0(a)  =  (r-l/r)/21/2;  (13) 

where  F  is  defined. as 

r  =  {(l+(l+2a2),/2l/21/2a}1/3  (13a) 

The  value  tQ  for  which  9*  =  0,  and  for  which  9'  reaches 
its  maximum  as  figure  3b  indicates,  is  a  decreasing 
function  of  a.  We  have  t0  =  »  for  a  -  0,  tg  *  l/2,/2 
for  a  =  l/23/2,  and  t0  =  0  for  a  =  <*>.  Furthermore, 
figure  4,  where  the  functions  t0(a),  t_(a)  and  t+(a)  are 
depicted,  shows  that  we  have  0<  t_  <  tg  <  t  +  <  “  for 
0  <  a  <  1/2372 


Fig.  4’—  The  Functions  t0(o),  t+(a)  and  t_(a) 


In  the  vicinity  of  the  points  of  stationary  phase  t± 
the  first  derivative  of  the  phase-function  may  be 
approximated  by  the  two-term  Taylor  series 
9’(t;a)  0?  (t-t±)©"  +  (t-t±)20"72,  (14) 

where  the  function  t±(a)  is  defined  by  equation  (12),  and 
9*and  9"'  represent  the  values  of  the  functions  9"(t;a) 
and  9"'(t;a)  for  t  =  t±(a).  By  using  equation  (12)  in 
equations  (lib)  and  (He)  we  may  obtain 
9’±=  ±2372a(I  -8»2)!/2/(l  +4a2±(l  -8a2)172]’72,  (15a) 
9”  =  T-217296a4/ll+4a2±(l-8e2)172J372.  (15b) 


In  the  neighborhood  of  the  point  tg,  where  9'=  0,  we 
have 

9'(t;e)  ~  0g  +  (t-t^'72,  (16) 

where  the  function  tg(a)  is  defined  by  equation  (13),  and 
9g  and  8g"  represent  the  values  of  the  functions  9'(t;a) 
and  9’"(t;a)  defined  by  equations  (Ua)  and  (He)  for  t  = 
tg(«).  The  functions  tg(a),t+(a)  and  t_(r),  9g(o)  and 
00 '(o),  0*>)  and  0'Ja),  8 '"(a)  and  0 ':'(»)  are  depicted 
in  figure  5.  The  signs  of  the  functions  0q,  8g',  8^.  and 
0JJ.'  are  readily  apparent  from  this  figure.  The  Taylor- 
series  approximation  (14)  and  (16)  are  useful  for  devising 
efficient  numerical  method  for  evaluating  the  integrals 
(1),  (3a, b)  and  (5a, b),  as  will  be  examined  in  detail 
elsewhere. 


Fig.  5  —  The  Functions  t0(a),  t±(o),  0g(a),  9g'(a), 

9^  (a)  and  0^'(a)  for  0<  a  <  0.5 

Asymptotic  approximations  valid  for  z  =  0  and  x 
-oo  will  now  be  obtained  for  the  integrals  (1),  (3a, b) 
and  (5a, b).  These  five  integrals  may  be  expressed  in  the 
form 

n+k(x,a)  =  Im^^+Cj^-),  (17) 

where  0  <  k  <  4,  Cj*  are  constants,  and  tpj*  are  the 
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fatam*  defined  m 

-  f"  exrti^^C^JKO^O)  dt;  (18) 

in  this  expression,  the  phase  0±(t;«)  is  defined  by 
•±(t;a)  -  (!7«tXl  +t2),/2  with  a  >0;  (19) 

finally,  the  constants  in  equation  (17)  and  the 
(tactions  ak(t)  in  equation  (18)  are  defined  in' the 


following  table: 

k  *k  h  «T 

0  ♦  1  11  (20a) 

1  ♦„  (l+t2)1/J  i  i  (20b) 

2  t(l+t2),/2  i  -i  (20c) 

3  ♦„  1+t2  -1  -1  (20d) 

4  *xy  t(l+t2)  -11  (20e) 


In  the  particular  case  a  =  0,  that  is  on  the  track  of 
the  ship,  equation  (18)  yields  ipk  =  ifk  =  yk-,  with 
Wk(x,°)  =  /0“  exp[ix(l+t2)1/2]  K(t)ak(t)  dt. 

The  major  contribution  to  this  integral  in  the  limit  x  -* 

-<*>  stems  from  the  point  of  stationary  phase  at  the  origin 
t=0.  Specifically,  we  may  obtain 
(n/2),/2(  -  x)1/2+k(x,0)  ~ 

Im  r^Kg  exp[i(x-it/4)]  as  x -* -<»,  (21) 

with  Cq  *  1,  Cj  =  i,  Cj  *  -1,  (21a, b.c) 

and  Cj  =  0  =  c4;  we  thus  have 
♦y  =  0  =  jxy  for  *  =  0,  (22a, b) 

in  accordance  with  the  symmetry  of  the  wave  pattern 
about  the  axis  a  =  0.  In  equation  (21),  Kg  represents  the 
value  of  the  function  K(t)  at  the  origin  t  =  0,  that  is  we 
have  Kg  =  K(0). 

For  0  <  a  <  l/23/2,  the  phase  0+(t;«)  is  stationary, 
that  is  B'+  =  0,  at  the  two  distinct  points  t_(o)  and  t+(e) 
defined  by  equation  (12),  whereas  0'_(t;o)  >  0  for  t  >  0, 
as  may  be  seen  from  figure  3b.  Equations  (17)  and  (18) 
then  yield 

n^k(x,cr)  ~  Im  as  x  -*  -o®  with  0  <  a  <  l/23/2. 

The  contribution  of  the  two  points  of  stationary  phase  t± 
can  be  evaluated  by  using  the  method  of  stationary  phase, 
with  the  result 

(«/2)l/2(-x)'%(x,a)  ~  Im  ck(Ak-E_  +Ak+E+) 
as  x  ■*  -°°,  with  0  <  a  <  l/23/2  and  (23) 

Cg  =  1,  c,  =  i  =  Cj,  c3  =  -1  =  c4;  (23a,b,c,d,e) 

furthermore,  Aj*  and  E±  are  the  amplitude  and 
exponential  functions  defined  as 


Ak±  =  K±a*/(T0"±)l/2, 

(24a) 

E±  =  exp[i(x0±  ±n/4)], 

(24b) 

where  K±,  aj*,  0±  and  9"±  are  defined  as 

K±=  K(t±),  ak*  =  ak(t±). 

(25a, b) 

e±  *  8(t±;«),  o;  «  e*<t±;.),  (2Jc,d) 

and  t±  is  the  function  of  «  given  by  equation  (12),  that  it 
we  have 

t±  =  [l±(l-8*2),/2)/4<».  (26a) 

We  may  then  obtain 

!+t±2  =  (1  +4«2±(I  -8«2)1/2]/8o2.  (26b) 

The  expressions  for  the  terms  ajf  may  readily  be  obtained 
from  equations  (20a-e),  (23b)  and  (26a, b).  Equations  (23c), 
(19)  and  (26a, b)  yield 

0*  =  (3  T(1  -8<r2)l/2]|l  +4^±(1  -8a2)l/2]1/2/2l/28«r.  (26c) 
Finally,  0’  is  given  by  equation  (13a),  that  is  we  have 
*0'±  =  23/2o(1  - 8o2)1/2/[1  +4o2±(1  -8o2)i/2]i/2.  (26d) 

It  may  be  shown  from  equations  (26a, c,d)  and 
verified  from  figures  4,  3a,  3c  and  3  that  we  have  t_  =  0 
and  0_  =  1  =  ©1  in  the  limit  o  -  0.  Furthermore, 
equations  (23b)  and  (20a-e)  show  that  we  have  ajT  =  1  for 
k  =  0,  1  and  3,  and  ak  =  O  for ).  =  2  and  4  in  the  limit 
a  =  0.  The  asymptotic  approximation  (23)  for  |k(x,o) 
therefore  becomes  identical  to  the  asymptotic 
approximation  given  by  equations  (21),  (21a, b,c)  and 
(22a, b)  for  $k(x,0)  in  the  limit  o  =  0,  if  the  contribution 
of  the  second  point  of  stationary  phase  t+  =  "is  null, 
that  is  if  Ak+  =  0  for  o  =  0  and  t  +  =  <*>.  In  other  words, 
the  asymptotic  approximation  (23)  for  |k(x,o)  is  uniformly 
valid  in  the  vicinity  of  the  track  of  the  ship  ar  =  0  if  Ak  = 
0  for  a  =  0.  Equations  (26a)  and  (26d)  yield  t  +  ~  l/2o 
and  -0^.  2a  as  a  -*  0.  We  then  have  -04  ~  l/t  + 
as  a  •*  0,  and  the  condition  for  the  asymptotic 
approximation  for  fk(x,o)  to  be  uniformly  valid  in  the 
limit  o  =  0  takes  the  form 

tl/2K(t)ak(t)  -  0  as  t  -  ~.  (27) 

In  the  limit  a  =  l/23/2,  we  have  '  =  l/2,/2  - 
t+  and  0^  =  0,  as  may  be  verified  from  equations 
(26a, d)  and  figure  5.  Equation  (24a)  then  shows  that 
have  |Aj*|  -*  00  as  a  -*  l/23/2;  and  the  asymptotic 
approximation  (23)  is  not  valid  in  the  vicinity  of  the 
boundary  of  the  Kelvin  wake.  A  complementary 
asymptotic  approximation,  expressed  in  terms  of  Airy 
functions,  valid  at  and  near  the  Kelvin  cusp  line  is  given 
in  Ursell  [11]  for  the  particular  case  of  a  pressure  point  at 
the  free  surface.  However,  we  are  mostly  interested  in  the 
sector  0  <  a  <  l/23/2,  that  is  inside  the  Kelvin  wake,  in 
the  present  study. 

The  far-fieid  asymptotic  approximation  (23)  shows 
that  the  wave  pattern  at  any  point  (x,a),  with  x  «  - 1 
and  0  <  a  <  l/23/2,  consists  in  two  elementary  plane 
progressive  waves.  Specifically,  equations  (4),  (20b),  (23) 
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and  (24b)  show  that  the  free-surface  elevation  within  the 
Kelvin  wake  at  a  suffidently-large  distance  behind  the  ship 
is  given  by 

(«/2)1/2(-x)l/2e(x,9)  ~  Re(A]exp(kpJ + A,+ expfxp + )), '  (28) 
where  the  phases  ip±  of  the  two  waves  are  gjvenby 
ip±  =x©±  ±it/4.  Curves  along  which  the  phases  or 
if  art  constant  then  are  defined  by  the  equation 

r 

x©±  ±n/4  =  constant.  This  relation  and  equation  (8)  then 
yield  the  following  parametric  equations  for  the  curves 
along  which  the  phase  is  equal  to  -2nn: 

-x±  =  (2nn±n/4)/©±,  (29a) 

y±  =  (2nn±n/4)«/e±,  (29b) 

where  0  <  a  <  l/23/2.  The  ten  constant-phase  curves 
corresponding  to  1  <  n  <  10  are  depicted  in  figure  6.  The 
“transverse”  and  “divergent"  waves  in  this  classical 
representation  of  the  Kelvin  wake  correspond  to  the  waves 
Afexp(i<p_)  and  A|+exp(kp+),  respectively,  in  equation 


Fig.  6  —  The  Classical  Kelvin  Ship  Wave  Pattern 


The  wavenumber  corresponding  to  the  wave  with 
phase  ip  ^  is  given  by  Vcp±.  The  corresponding 
wavelength,  say  A±,  and  direction  of  propagation  with 
respect  to  the  track  of  the  ship,  say  0±,  then  are  given  by 
A±  =  2n/|V<p±|  and  p±  =  tan-l(<p*/(p*),  where  <p*  and 
cp*  represent  the  x-  and  y-derivatives  of  cp±  and  |V(p±| 

=  l(<p* )2  +  (<P* )21,/2.  The  relation  <fx  -  x9±  ±n/4  and 
equation  (26c)  then  yield 

A±  =21/216no2/[3?(l-8o2)1/2Hl-4o2±(l-8«2)1/:]1/2,  (30a) 
P±  =  sitf1{(l-4n2±(l-8«2)1/2l/[l+4<»2±(l-8o2)w2i}1/2.  (30b) 
Equations  (30a, b)  show  that  we  have  A_  =  2n,  p_  ~  0 
and  A+  =  0,  /3+  =  n/2  in  the  limit  a  =  0,  and  k_  - 
4n/3  =  A+>  /?_  =  sin- l(l/3,/2)  =  0+  for  a  =  l/23/2. 
More  precisely,  we  have 

2n>l_>  4n/3  >  A+  >  0  and  (31a) 

0  <  p_  <  sin"  ‘(l/3,/2)  <  p+  <  n/2,  (31b) 

as  may  be  seen  from  figure  7  where  the  functions 
A±(a)/2n  and  2p±(a)/n  are  depicted. 


Fig.  7  —  The  Wavelengths  A±(o)  and  Propagation  Angles 
P±(a)  of  the  Transverse  and  Divergent  Waves  in  the 
Kelvin  Wake 


Equation  (28)  shows  that  the  amplitudes  of  the 
transverse  and  divergent  waves  in  the  Kelvin  wake  are 
asymptotically  given  by  (2/n)1/2|Ajt|/(-x)1/2  as  x  -* 

The  steepnesses,  say  s±,  of  these  waves  then  are  given  by 
g±»  (2/n)1/2|A|±  |/(  -  x)i/2A  ±  •  Equations  (24a),  (25a, b) 
and  (20b)  then  yield 

(-x),/2s±  ~  o±(o)|K(t±)|  as  x  -  (32) 

where  o±(o)  is  defined  as  (2/n)1/2(l+t2.)1/2/(:F0±)1/2A;|;. 

Equations  (26b, d)  and  (30a)  then  yield 
u±  =  (3t(1  -8o2)1/2](1  -4o2±(l  -8o2),/2J1/2 
[1  +4a2±(l  — 8a2)*/2]3/V 

64n(2,/2n)1/V/2(l-8o2)IM.  (33) 

Equations  (12)  and  (33)  yield  t_  =  0  and  o_  = 
l/n(2n)1/z  in  the  limit  a  =  0,  for  which  we  have  A_  =  2n 
and  p_=  0  as  was  noted  previously.  The  steepness  of  the 
transverse  wave  at  a  point  (x,0)  on  the  track  of  the  ship 
then  is  given  by 

s  _  (x,0)  ~  |K(0)|/n(2n)l/2(  -  x),/2  as  x  -  -oo.  (34) 

Equations  (12)  and  (33)  also  yield  t+  ~  l/2o  and 
o+<v  l/16n3/2o7/2  in  the  limit  o  -*  0,  for  which  we  have 
A+-*  0  and  P+~*  n/2.  We  thus  have  o+~  t7+/2/n(2n)1/2  as 
a  **  0,  and  equation  (33)  shows  that  the  steepness  of  the 
divergent  wave  at  a  point  (x,o)  in  the  vicinity  of  the  track 
of  the  ship  is  given  by 
s+(x,o)  «v  t7+/2|K(t+)|/n(2n)1/2(-x)1/2 
as  a  —  0,  with  t+  ~  l/2o.  (35) 

The  steepness  s+(x,a)  then  becomes  unbounded  as 
a-Oif 

t7/2|K(t)|  -*  *  as  t  -*  oo.  (36) 
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Equations  (4),  (20b)  and  (27)  show  that  the  asymptotic 
expansion  (23)  for  the  free-surface  elevation  is  valid  in  the 
vicinity  of  the  track  of  the  ship  if 
t3/2K(t)  -  0  as  t  -  co.  (37) 

Let  us  assume  that  we  have 

|K(t)|  -v  l/tH  as  t  -  oo.  (38) 

Both  conditions  (37)  and  (36)  are  then  satisfied  if 
3/2  <  ft  <  7/2.  (39) 

Condition  (37)  means  that  the  amplitudes  of  the 
divergent  waves  in  the  Kelvin  wake  vanish  as  a  —  0,  that 
is  as  the  track  of  the  ship  is  approached,  whereas 
condition  (36)  means  that  the  waves  become  infinitely 
steep;  this  is' theoretically  possible  because  the  wavelengths 
of  the  divergent  waves  vanish  as  a  -*.0,  as  is  indicated  in 
equation  (31a)  and  figure  7.  However,  infinitely-steep 
water  waves  cannot  exist  in  reality;  indeed,  there  exists  a 
theoretical  upper  bound  for  the  steepness' of  water  waves 
in  deep  water  which  is  approximately  equal  to  1/7. 
Condition  (39)  therefore  suggests  that  noy  divergent  waves 
can  exist  within  a  certain  domain  in  the  vicinity  of  the 
track  of  the  ship,  and  that  the  Kelvin  wake  contains  three 
distinct  regions:  (i)  an  inner  region  adjacent  to  the  track 
of  the  ship  where  only  transverse  waves  can  exist  as. was 
just  noted,  (ii)  an  outer  region  where  both  transverse  and 
divergent  waves  are  present,  and  (iii)  a  region  at  the 
boundary  between  the  inner  and  outer  regions  where  short 
steep  divergent  waves,  as  well  as  transverse  waves,  can  be 
found.  It  must  be  kept  in  mind,  however,  that  these 
conclusions  regarding  the  Kelvin  wake  are  based  on 
condition  (39),  which  was  obtained  on  the  basis  of  an 
analysis  in  which  surface-tension  and  nonlinearities  are 
ignored.  Inasmuch  as  this  linear  no-surface-tension 
analysis  predicts  short  steep  waves,  both  surface-tension 
arid  nonlinear  effects  are  liable  to  be  significant,  and  these 
effects  should  therefore  be  included  in  a  more  realistic 
analysis.  In  particular,  it  is  evident  from  Lamb  (9, 
pp.  468-470)  and  Wehausen  and  Laitone  (10,  pp.  636-637) 
that  the  system  of  divergent  waves  in  the  immediate 
vicinity  of  the  track  of  the  ship  may  be  profoundly 
affected  by  surface  tension. 


Fig.  8  •—  Waterlines  and  Framelines  of  the  Simple  Ship 
Bow  Form  Considered  for  Numerical  Applications 


A  SIMPLE  TEST  CASE: 

THE  FAR-FIELD  WAVE-AMPLITUDE  FUNCTION 
The  foregoing  theoretical  results  are  investigated 
numerically  for  the  simple  semi-infinite  ship  form  studied 
previously  by  Scragg  (5).  This  ship  form  consists  of  a  bow 
region,  with  length  L,  followed  by  a  parallel  body,  with 
invariant  framelines,  extending  to  infinity  downstream.  All 
framelines,  both  in  th:  bow  region  and  downstream  from 
it,  are  trapezoidal  in  shape  with  constant  draft  D.  The 
waterlines  are  rectilinear.  More  precisely,  the  hull  form  is 
defined  by  the  equations 
y-  ±  (tan/3 +z  tany)(l  -x) 

for  0  <  x  <  1  and  0  >  z  >  -d,  (40a) 

y  =  ±(tan/3  +  z  tany)  for  x  <  0  and  0  >  z  >  -d,  (40b) 

where  x,y,<vand  d  are  nondimensional  in  terms  of  the 
length  of  the  bow  region,  that  is  we  have7=X/L  and 
d  =  D/L.  Equations  (40a, b)  require  that  the  condition 
tan/3  >  d  tany  (41) 

be  satisfied.  Equation  (40a)  shows  that  the  entrance  angle 
at  the  bow  (x  =  1,  z  =  0)  is  equal  to  2/3,  and  it  may  be 
seen  from  equations  (40a, b)  that  y.  represents  the  flare 
angle  for  x  <  0.  The  four  waierlines  corresponding  to 
z  =  0,  -d/3,  -2d/3,  -d  and  the  five  framelines 
corresponding  to  x  =  1,  0.75,  0.5,  0.25,  0  are  depicted  in 
figure  8  for  d  =  0.1,  p  -  12°  and  y  =  45°.  The 
notation 
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Po  =  ta np,  r0  =  lany  <42a,b) 

will  be  used  for  shortness  hereafter. 

The  far-field  wave-amj  litude  function  K(t)  for  the 

foregoing  semi-infinite  hull  has  been  evaluated  for  two 

simple  approximations  defined  explicitly  in  terms  of  the.. 

hull  shape  and  the  Froude  number,  namely  the  Michell , 

'  » 

thin-ship  approximation  and  the  zeroth-order  slender-ship, 
approximation  (7).  for  which  the  function  K(t)  is  denoted 
KM(t)  and  K0(t),  respectively.  The  Michell  thin-ship 
approximation  is  given  by  the  product  of  two  jingle 
integrals,  as  follows: 

Kn,  =  -v*  /q  cxp[-i«r(l  +  t2)l/2x]  dx 

/d°  exp[^(l  +t2)z]  (Po+YqZ)  dz. 

t. 

These  integrals  can  be  evaluated  analytically,  with  the 
result 

KN,(l)  =  4I/?0 -  f/J0 - r0d)e -  F2y0(l  - e)/(l  + 12)} 

sin[v2(l  +  t2)l/2/2|  expl-i^O  +t2),/2/2]/(l  +t2)3/2,  (43) 
where  the  term  c  is  defined  as 

c  =  cxpi-v2d(l  +t2)].  (44) 

Equations  (43)  and  (44)  yield 

|KM(t)|  ~  4/l0|sin(v2t/2)l/t3  as  t  -  °°.  (45) 

The  zeroth-order  slender-ship  approximation  may  be 
expressed  in  the  form 

K0(t)  =  Kq  (t)  +  Kq  (t),  (46) 

where  K^(t)  is  given  by  the  suin  of  a  double  integral  over 
the  hull  surface  and  a  single  integral  along  the  top 
waterline,  as  follows: 

K*  =  vA/0l  dx  f  ^  dz  cxp[v2(l  +t2)z] 
cxpl  —  iv2(l  +  t2),/2(x  ±  yt)l(/30 + y0z) 

~v~  dx  cxpl  —  iv2(  1  +t2)1/2(x±yt)] 

/33/[l  +/?5  +  y2(l  -x)2J.  (47) 

The  integration  in  these  integrals  is  carried  over  the 
positive  half  of  the  hull  surface.  Equations  (40a)  and 
(42a, b)  then  yield  y  =  (/J0  +  y0z)(l  -x).  We  then  have 
x  +  yt  =  I  -  (1  i/)0t)(l  -x)±y0t(l  -x)z. 

By  using  this  relation  into  equation  (47)  we  may  then 
obtain 

cxpliv2(l+t2)'  2|K0’(t)  =  v2Jq  cxp|iv2(l  +t2),/2(l  T/J0t)x) 
l/Vo  ■  y„l f  -P^(l+^+y5x2))  dx,  (48) 


where  I*,  with  n.=  0  and  1,  are  the  integrals  defined  as 

(-  Iff*  =  v2//  apl^d  +t2Xl  +iy0ux)z]zn  dz, 
with  u  defined  as 

u  =  t/(l+t2),/2.  (49) 

The  integrals’ll  may  be  evaluated  analytically  with 
the  result 

I =  (1  -«Cj.)/(l  +t2Xl  +iy0ux). 

If  =  [F2(!-ee£)/(I+t2Xl?iy0ux)-deeJ 
/(I  +t2XI  Tiy0ux), 

where  e  is  the  exponential  function  given  by  equation  (44) 
and.e.,  is  the  exponential  function  defined  as 
e.,  =  expl±iv2dy0t(l+t2)1/2x).  (50) 

We  then  have 

Po't?  ~Yol?  +/Jo+)'ox2)  = 

At  -eesB±/(l  +t2Xl  +  iy0ux),  (51) 

where  the  terms  A„  and  B,  arc  given  by 
A±  =  /i0/(l  +t2Kl  +  ir0ux)-/ij*/(l  +Pq+ TqX2) 
-y0F2/(l+t2)2(lTiy0ux)2, 

B  ±  =  Po~Yod~  r0F2/0  +  »2Xi  *  >y0ux). 

It  may  be  verified  that  we  have 

A,  =  P0l(l-^t2)/(l+l2)  +  x2/)2y2/(I+/i2  +  y2))/(l+/)2) 

+  r0C±,  (52a) 

Bt/(l+t2XlTiy0ux)  =  /ya+t2)-/,)0;..  (52b) 

where  the  terms  and  D±  arc  given  by 

c±  =  IP2rJ/(i+^Xi+^+yg)]x2d-x2)/(i+/)g+rgx2) 

-  l/Vo'2*2 + F2(  1 + 12  -  yo*2x2)/(i + 12 + y^x2) 

Tit(l  +t2),/2{/)0-2FV(I  +t2  +  y2t2x2)}x] 

/( 1  + 12)(  1  + 12 + ygt2x2),  (53a) 

(1  + 12)(  1  + 12  +  yJl2x2)D  ±  =  /)0y0t2x2  +  d(l+t2) 

+  F2(1  +t2-yQt2x2)/(l  + 12  +  y^t2x2) 

T  it(l  +  t2)1,2I/i0-  y0d  -  2F2y0/(l  + 12+ y*t2x2)]x.  (53b) 
Equation  (50)  yields 
cxp[iv2(l  +t2)l/2(l+/)0t)x]c1.  = 

cxp[iv2(l  +t2),/2(l  T/)dt)x],  (55) 

where  pd  is  defined  as 

Pd  =  /)0-dy  o  =  lanJJ-dtany,  (56) 

as  may  be  obtained  from  equations  (42a, b).  By  using 
equations  (55),  (51),  (52a, b)  and  (44)  into  equation  (48)  we 
may  obtain 

cxp[iv2(l  +t2)l/2)K(f  (t)  = 

PotO -/}o,2)(|+t2)  ''o-O^+^oC  +Po+r5)-'lf«;P0) 

-  (I  +  /$cxp{  -  v2d(l  +  t2)}(l  +t2)-  'I* (t;/)d)]/(l  +/)2) 

+  y0(J  *  C;Po-c  ±)  +  cxp{  -  v2d(l  + 12)}  J  4  (t  ;pd,D  ± )),  (57) 

where  the  functions  I*(t;/i)  and  J^Ot/i.A),  or  more 
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precisely  I* (t #.v)  and  J±(t;/).A,v),  are  defined  in  terns of 
the  integrals 

I* m  =  v2^1  exp^l  + t2),/2(l  TjOxJx23'  dx.  (51a) 

J±(t^,A)  =  v2/,1  exppvty  +t2),/2(l  7Jh)x]A(x)  dx.  (51b) 

The  integrals  I*  0$)  can  be  evaluated  anslyticaly, 
with  the  result 


lfM)  =  i(l-E±)/o, 


(59a) 


lf(»;fl)  =  -  i((l + 2iF2/o±  )£±+2F4(l-Etya2t]/ot.  (59b) 
where  and  E,  are  defined  as 
ot  =  (l+t2)uHl?pl).  E±  *  cxp0^o£).  (59c.d) 

Expressions  (59a,b)  for  ig  (t'JI)  and  l,+(t^)  art  not  valid 
in  the  special  case  when  we  have  t  <=  l/p,  for  which 
expressions  (59a,  b)  become 

Iq"  (i//hp)  =  v2,  i;a/m  =  Jn.  (59e,n 

The  integrals  J±  defined  by  equation  (58b),  where 
the  amplitude  function  A(x)  lakes  the  form  of  C±  or  D± 
specified  in  equations  (53a,b),  cannot  be  evaluated 
analytically.  These  integrals  were  then  evaluated 
numerically  by  c’viding  the  integration  range  0  <  x  <  I 
into  N  segments  of  equal  length  and  using  piecewise 
quadratic  approximations  for  the  amplitude  function  A(x) 
within  each  segment.  In  this  manner,  we  may  obtain.  . 

N 

o.JJttp.A)  ~Z(£i)'-,ii(Aj-t±AJtl) 

+  (NF2/o  ±  Kt  t  -  1  KAj  H  -  Aj) + (2NF2/o  £ ) 

{ 1  + 1 1  +  it2NF2/o  ,K£t-  1)}(Aj  +Aj+|-2Aj+  |/2)J,  (60) 
where  o±  is  given  by  equation  (59c),  is  defined  as 
i±  =  exp(iv2oI/N),  (60a) 

and  Af,  AJ  + ,,  AJ+  ,/2  represent  the  values  of  the 
amplitude  function  A(x)  at  the  points  Xj  =  (j-l)/N,  Xj+, 
=  j/N,  and  xf ,  |/2  =  (j-l/2)/N,  respectively.  Expression 
(60)  for  J  t(t;/5,A)  is  not  valid  in  the  special  case  when  I 

=  1/p,  for  which  wc  have 
N 

6NJ  ( (I //);/), A)  AJ  +  AJt,+4AJ  +  1/2.  (60b) 

In  summary,  the  zeroth-order  slender-ship 
approximation  K0(t)  is  determined  by  equations  (46),  (57), 
(59a-f),  (60),  (60a, b),  (42a,b)  and  (56).  In  the  limiting  case 
y0  =  0,  these  equations  yield 
K0(()  =  4/3o(l+t2r3'2((l-/$2)(l+02)-' 

-  exp{  -  v2d(l  + 12)} ](R  -  il)/(l  -flgt2),  (61) 

where  the  terms  R  and  I  are  given  by 
R  ■=  sipjv2(l  +  l2)l/2/2)  coslyfy  -M2)1/2'& 

P, 


la  the  thin-ship  limit  p0  —  0,  equations  (61)  and  (61a,b) 
yieM 

K</t)  ~  4^1  +t2)-3/2(l -exp{  -i^dd  +t2)}]. 

stnfv^l  +t2),/2/2]  expf-i^l  +t2)1/2/2J  as  p0  -  0.  (62) 

The  thin-ship  limit  (62)  of  the  zeroth-order  slender-ship 

approximation  Kq(0  may  be  seen  to  be  identical  to  the 

Michel  thin-ship  approximation  (43)  in  the  particular  case 

f0  *=  0.  However,  the  thin-ship  limit  (62)  of  the  slcnder- 

*P  approximation  (61)  is  not  uniformly  valid  in  the  limit 

t  -*  •»:  indeed,  equations  (61)  and  (61a,b)  yield  equation 

(62)  in  the  limit  p0 1  -*■  0.  Equation  (62)  yields 

ItytH  ~  4/J0|sin(v2t/2)|/t3  for  1  «  t  «  1%.  (63) 

More  generally,  equation  (61)  yields 

lKo(t)|  *v  4^0|R-iI|/(l  +/J2)!3  as  t  -  ».  (64) 

Equations  (61a, b)  show  that  |R-il|  ~  /)0t|sin(v2/)0t2)|/2  if 

Pfj,  »  1.  We  then  have 

IK^OI  ~  2^|sin(v2/30t2)i/(l  +pftl2  for  t  »  !//)„.  (65) 

In  the  liptit  t-*»,  equations  (57)  and  (59a-d)  yield 
'  exptfv^O)  ~  /)0[(i  -pfr2)l0*(i;p0) 

+rfrl  +Pl+ry  ,/7gt2I,i(t;fJ&)]/(l  +  p20) t2  as  t  - 
where  we  have 

tfOlW  ~  i[!-exp{iv2t(lT/30t)}l/(l  +  /)0t)t  as  t  - 
IfOlPo)  ~  -*  expjivfyl  ?/J0t)}/(l  */)0t)t  as  t  -  <=°. 

We  then  have 

fiyOl  ~  2/J0|N|/(l  +  /)2)| I  - /)2i2|t3  as  t  -  oo,  (66) 

where  the  term  N  is  given  by 
N  =  (1  —  /Jot2)Isin(w2i)  +  i  cos(v2t)) 

-0  -°2Poi2)\P(f  sin(v2/30t2)  + 1  cos(v2/J0t2)],  (66a) 

with  o2  defined  as 
,2  -  ri  a.ni\/n  j. „2\ 


«*  =  d+/^y(i  +Po+y& 


(66b) 


We  may  then  obtain 
'  |Ko(0l  ~  4/30|sin(v2t/2)|/(l+/32)l3  for 
I  «  t  «  F/(Pq),/2,  (67a) 

|K0(t)|  ~  2/$sin(v2f?0t2)|/(l+/32  +  y2)t2  for  (  »  1 %.  (67b) 
Equations  (67a, b)  arc  identical  lo  equations  (63)  and  (65) 
in  the  limits  p0  -*  0  and  y0  -*  0,  respectively.  More 
generally,  equation  (66a)  yields  the  following  upper  bound 
for  the  term  |N|  in  equation  (66) 

INI2  <  (l-/$2)2  +  (l+/$2)(l  O2/)2!2)2 


+  2(1  +  /)0t)2(l  +o/?0l)|(1  -  /30t)(  1  o/)„t)|. 


(68) 


Equations  (67a, b)  thus  show  that  in  the  limn  i  —  «> 
we  have  |K0(t)|  ~  1/t3  for  I  «  t  «  F/(/)0)'  2,  that  is  lor 
moderately  large  values  of  t,  and  th0(Dl  ~  Ur  for  I  pn 


.<&  t,  th(il  is  for  very  large  values  of  i  The  asymptotic 
+  «Pjv*(l  +  P)'  ‘/2J  GKlv%t  (b/<r ,b)  and  (45)  show  that  the  Michcll 

'o'  +  +<fy  %  .  *»,  -  KM(t)  corresponds  to  the  thin-ship  limn  pn 

I  *  sin  (v^l  +  t2)*  2/2| - iinfyyfyjlO  + l2}1  2/J).  «  )  of  the  zeroth-order  slender-ship  approximation  K0(t), 
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Copy  available  to  DTIC  does  not 
permit  fully  legible  reproduction 


Fig.  9a  —  Real  and  Imaginary  Parts  of  the  Functions 
KM(t)  and  K0(t)  for  a  Simple  Ship  Bow  Form  with 
p  =  12°  (Top)  and  3°  (Bottom),  y  -  0  and  F  =  0.5 


and  that  this  thin-ship  limit  is  not  uniformly  valid  in  the 
limit  t  -*  oo.  More  generally,  the  limiting  processes  p0  •*  0 
and  t  -*  «  cannot  be  interchanged. 

Figures  9a,b  depict  the  real  and/or  imaginary  parts 
of  the  functions  (1  +  t2)’  2KM(t)  and  (1  +  t2)|/2K0(t)  for  the 
simple  bow  shape  defined  by  equations  (40a, b),  where  the 
nondimensional  draft  d  and  the  maximum  flare  angle  y 
arc  taken  equal  to  d  =  0.1  and  y  =  0,  and  four  values 
of  the  half-cntrancc-angle  ft  arc  considered,  namely  ft  = 
12°,  3°,  P  and  20'.  The  Froude  number  based  on  the 
length  of  the  bow  region  is  taken  equal  to  F  =  0.5  in  the 
numerical  results  presented  in  figures  9a,b  and  in  figures 
10  and  1 1  considered  further  on.  The  values  of  l//30  = 
1/tan/)  corresponding  to  the  values  of  ft  equal  to  12°,  3®, 
P  and  20  are  approximately  equal  to  4.7,  19,  57  and 
172,  respectively.  The  functions  (1  +t2),,2KM(t)  and 
(1  +t2)'  2K0(t)  arc  depicted  for  0  <  t  <  lo  in  figures  9a, b. 
Differences  between  the  approximations  KM  and  K0, 
especia'ly  their  imaginary  parts,  can  be  seen  to  be 
s  ibstaitial  in  figure  9a  corresponding  to  p-  12°  and  3°. 
Figure  9b  shows  that  differences  between  the  imaginary 
parts  of  the  functions  KM  and  K0  remain  appreciable  even 


Fig.  9b  —  Real  and  Imaginary  Parts  of  the  Functions 
KM(t)  and  Ko(t)  for  a  Simple  Ship  Bow  Form  with 
p  =  1°  (Top)  and  20’  (Bottom),  y  =  0  and  F  =  0.5 


for  values  of  ft  equal  to  1°  and  20',  which  arc  quite 
small,  and  for  values  of  t  that  arc  much  smaller  than 
1/P0. 

The  top  and  bottom  parts  of  figure  10  depict  the 
real  and  imaginary  parts  of  the  functions  (1  +  t2)|/2KM(t) 
and  (1  +t2)|/2K0(t),  respectively,  for  the  previously- 
considered  simple  bow  shape  with  d  =  0.1,  p  =  1 2°,  F 
=  0.5  and  for  two  values  of  the  flare  angle,  namely  fo;  y 
=  0  and  45°.  The  top  part  of  the  figure  shows  that 
differences  between  the  curves  corresponding  to  y  =  0 
and  y  -  45°  arc  faily  small,  and  are  appreciable  only  for 
small  values  of  t,  for  the  Michell  thin-ship  approximation 
Km.  In  particular,  the  asymptotic  approximation 
KM(t)  ~  4p0  sin(v2t/2)  exp(-iv2l/2)/t3  as  t  —  00 , 
which  may  be  obtained  from  equations  (43)  and  (44),  is 
independent  of  y.  The  bottom  part  of  figure  10  shows 
differences  between  the  curves  corresponding  to  y  =  0 
and  y  -  45°  for  the  slender-ship  approximation  K0  that 
are  significantly  larger  than  those  for  the  Michell 
approximation  KM,  especially  for  intermediate  values  of  t 
in  the  vicinity  of  t  -  l/tan/)  ~  4.7.  The  flare  angle  y 
thus  has  a  pronounced  effect  upon  the  behavior  of  the 
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0.03 


Fig.  10  —  Real  and  Imaginary  Parts  of  the  Functions 
KM(t)  and  K0(t)  for  a  Simple  Ship  Bow  Form  with 
y  =  0  and  45°,  ft  =  12°  and  F  =  0.5 


slender-ship  approximation  KqO)  for  values  of  t  in  the 
neighborhood  of  l/tan/3.  The  influence  of  y  upon  |Kg(t)| 
for  large  values  of  t  is  explicitly  indicated  by  the 
asymptotic  approximations  (67a, b). 

Finally,  figure  1 1  depicts  the  real  and  imaginary 
parts  of  the  functions  (1  +  t2),/2KM(t)  and  (1  +  t2)1/2K0(t) 
for  the  simple  bow  shape  depicted  in  figure  8,  for  which 
we  have  d  =  0.1,  p  =  12°  and  y  =  45°,  at  a  value  of 
the  Ffoude  number  F  equal  to  0.5.  Differences  between 
the  approximations  KM  and  K0  may  be  seen  to  be  quite 
substantial.  In  particular,  the  function  (1  +  t2)|/2K0(t)  has 
a  peak  at  t  =  1 /tan/1  =  1  /pQ. 


p  »  12*.  -y  «  45*.  F  »  0.5 


Fig.  1 1  —  Real  (Top)  and  Imaginary  (Bottom)  Parts  of 
the  Functions  KM(t)  and  Ko(t)  for  a  Simple  Ship  Bow 

Form  with  p  =  12°,  y  =  45°  and  F  =  0.5 

A  SIMPLE  TEST  CASE: 

THE  FAR-FIELD  KELVIN  WAKE 

The  expressions  for  the  far-field  wave-amplitude 
function  K(t)  obtained  in  the  foregoing  section  for  a 
simple  bow  shape  and  for  two  simple  approximations  to 
the  function  K(t),  namely  the  Michell  thin-ship 
approximation  KM  and  the  zeroth-order  slender-ship 
approximation  K„,  may  now  be  used  into  the  previously 
determined  asymptotic  approximations  for  the  far-field 
Kelvin  wake. 

Far  behind  the  ship,  that  is  for  x  —  -“>,  the 
functions  jk(x,o),  where  0  <  k  <  4,  defined  by  equations 
(17),  (18),  (19)  and  (20a-e)  are  given  by  the  asymptotic 
approximation  (23).  The  real  and  imaginary  parts  of  the 
amplitude  functions  Ak  (a)  and  A^  (a)  in  this  asymptotic 
approximation,  given  by  equation  (24a),  arc  depicted  in 
figures  12a  and  b,  13a  and  b  and  14  for  0  <  a  <  l/2}/2. 
More  precisely,  figures  12a, b  and  13a, b  represent  the 
amplitude  functions  Ay  and  Ak+  for  k  =  0,1, 2, 3 
associated  with  the  potential  $  and  its  derivatives  4X,  iy 
and  4XX.  and  correspond  to  the  approximations  KM  and 
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Kg,  respectively.  Figure  14  depicts  the  amplitude  functions 
A4"  and  A4  corresponding  to  4xy  for  the  approximations 
Km  and  Kg. 


Fig.  12a  —  Real  (Top)  and  Imaginary  (Bottom)  Parts  pf 
the  Functions  Kta*A*®jt)l/2  for  k  =  0  (Left)  and  1 
(Right),  K  =  Km  (Michell  Thin-Ship  Approximation)  and 
a  Simple  Ship  Bow  Form  with  p  =  12°,  y  =  45°  and 
F  =  0.5 


Fig.  12b  —  Real  (Top)  and  Imaginary  (Bottom)  Parts  of 
the  Functions  K  taf /(TB'J)^2  for  k  =  2  (Left)  and  3 
(Right),  K  =  Km  (Michell  Tlnn-Ship  Approximation)  and 
a  Simple  Ship  Bow  Form  with  p  =  12°,  y  -  45°  and 
F  =  0.5 

It  may  be  seen  from  figures  12a  and  b,  13a  and  b 
and  14  that  the  amplitude  functions  Ak~  (cr)  and  Ak  (a) 
become  unbounded  in  the  limit  a  -*  1/23'2.  This 
singularity  at  o  =  1/23  '2  stems  from  the  fact  that  the 
asymptotic  approximation  (23)  is  not  uniformly  valid  in 
the  limit  a  —  1/23'2,  as  was  already  noted.  The  amplitude 
functions  Ak+  to)  corresponding  to  the  system  of  divergent 
waves  in  the  Kelvin  wake  also  become  unbounded  in  the 
limit  a  -»  0  in  figures  13b  and  14  corresponding  to  k  =  2 


Fig.  13a  —  Real  (Top)  and  Imaginary  (Bottom)  Parts  of 
the  Functions  K^/fTe")^2  for  k  =  0  (Left)  and  1 
(Right),  K  =  Kg  (Zeroth-Order  Slender-Ship 
Approximation)  and  a  Simple  Ship  Bow  Form  with 
p  =  12°,  y  =  45°  and  F  =  0.5 


the  Functions  K±af /(T  9±),/2  for  k  =  2  (Left)  and  3 
(Right),  K  =  Kg  (Zeroth-Order  Slender-Ship 
Approximation)  and  a  Simple  Ship  Bow  Form  with 
P  -  12°,  y  =  45°  and  F  =  0.5 

and  3  for  the  slender-ship  approximation  K0  and  to  k  =  4 
for  both  approximations  KM  and  K0,  respectively  This 
singularity  at  a  =  0  illustrates  the  previously-noted 
conclusion  that  the  asymptotic  approximation  (23)  is  not 
uniformly  valid  in  the  vicinity  of  the  track  of  the  ship  a 
=  0  if  condition  (27)  is  not  satisfied.  Equations  (20a-c) 
show  that  we  have  ag  ~  1,  a,  ~  t,  a2  ~  t2,  a3  ~  t2  and 
a4  ~  t3  as  t  -*  oo,  and  equations  (45)  and  (67b)  yield 
|Km|  ~  1/t3  and  |Kg|  ~  1/t2  as  t  -*  »,  respectively. 
Condition  (27)  therefore  is  not  satisfied  for  k  >  4  and  k 
Sr  2  for  the  approximations  KM  and  Kg,  respectively 
Condition  (27)  however  is  satisfied  for  k  =  0  and  1, 
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corresponding  to  the  potential  |  and  the  free-surface 
elevation  |x,  for  both  the  approximations  KM  and  Kq. 

Figures  12a  and  13a  show  that  the  amplitude 
functions  and  Af,  corresponding  to  the  transverse 
waves  in  the  Kelvin  wake,  generally  are  larger  in 
magnitude  than  the  amplitude  functions  A^  and  A,"*" 
corresponding  to  the  divergent  waves,  whereas  the  reverse 
may  generally  be  seen  to  hold  in  figures  12b  and  13b  for 
the  amplitude  functions  Af  and  A  * .  The  relative 
importance  of  the  divergent  waves  with  respect  to  the 
transverse  waves  thus  increases  with  k.  Indeed,  the 
transverse-wave  amplitude  function  A f(a)  is  hardly  visible 
on  the  scale  of  the  divergent-wave-amplitude  function 
Af  (a)  used  in  figure  14. 


a  - ir,  t  -  ar,  r  -  oj 


Fig.  14  —  Real  (Top)  and  Imaginary  (Bottom)  Parts  of 
the  Functions  K±a.f7(  T9X)1/2  for  KM  (Left)  and  K0 
(Right)  and  a  Simple  Ship  Bow  Form  with  p  =  12°, 
y  =  45°  and  F  =  0.5 

The  divergent-wave-amplitude  functions  Ak+  (a) 
associated  with  the  approximation  K0  are  most  notably 
different  from  the  corresponding  functions  Ak+(o) 
associated  with  the  Michell  approximation  KM  for  values 
of  a  in  the  vicinity  of  a  =  0,  as  was  already  noted,  and 
of  a  =  tan/)/(2  +  tan2/)).  In  the  vicinity  of  this  value  of  a, 
the  divergent-wave-amplitude  functions  Ak+  (a)  associated 
with  the  approximation  Kq  exhibit  a  peak,  which  is  quite 
pronounced  for  k  >  1.  The  foregoing  particular  value  of 
a  corresponds  to  the  special  case  when  the  point  of 
stationary  phase  t  +  ,  defined  by  equation  (12),  is  equal  to 
the  value  l/p0  =  1/tan/)  for  which  the  function  Kfl(t) 
displays  a  peak,  as  may  be  seen  from  figures  10  and  II. 

Figure  15  depicts  the  amplitude  functions 
(2/u)l/2A|_  and  (2/n)l/2A,+  that  are  associated  with  the 
free-surface  elevation  fai  behind  the  ship,  as  is  specifically 


2  +  tan** 


Fig.  15  —  Amplitude  of  the  Transverse  (-)  and 
Divergent  ( ^  )  Waves  in  the  Kelvin  Wake  for  a  Simple 
Ship  Bow  Form,  with  /)  =  12°,  y  =  45°  and  F  =  0.5,  in 
the  Michell  Thin-Ship  Approximation  (Left)  and  the 
Zeroth-Order  Slender-Ship  Approximate  i  (Right) 

indicated  in  equation  (28),  for  the  previously-considered 
ship  bow  shape.  It  may  be  seen  that  the  amplitude  a  +  (o) 
of  the  divergent  waves  in  the  Kelvin  wake  vanishes  as  a  -* 
0,  that  is  at  the  track  of  the  ship,  and  is  generally  smaller 
than  the  amplitude  a_(o)  of  the  transverse  waves;  this  is 
especially  true  in  the  vicinity  of  the  track  of  the  ship. 
Differences  between  the  wave-amplitude  functions  a±(<r) 
corresponding  to  the  approximations  KM(t)  and  Xc(t)  are 
particularly  striking  for  the  amplitude  a+(u)  of  the 
divergent  waves  in  the  vicinity  of  a  =  tan/)/(2+tan2/)), 
where  the  function  a+(o)  associated  with  the 
approximation  K^t)  exhibits  a  sharp  peak. 

Figure  16  depicts  the  wave-steepness  functions  s_(o) 
and  s+(o),  which  correspond  to  the  ratios  of  the  wave- 
amplitude  functions  a _(<r)  and  a  +  («)  depicted  in  figure  15 
over  the  wavelength  functions  A  _  (a)  and  A  +  (o)  defined  by 
equation  (30a);  the  steepness  functions  s±(o)  are 
specifically  defined  by  equations  (32)  and  (33).  The 
divergent  waves  in  the  Kelvin  wake  may  be  seen  to  be 
generally  steeper  than  the  transverse  waves,  even  though 
figure  15  shows  the  transverse  waves  to  be  larger  in 
amplitude  than  the  divergent  waves.  This  is  especially  true 
near  the  track  of  the  ship  where  the  steepness  of  the 
divergent  waves  becomes  infinitely  large,  even  though 
figure  15  shows  that  their  amplitude  vanishes  as  o  -*  0. 
The  divergent  waves  become  infinitely  steep  at  the  track 
of  the  ship  because  the  wavelength  A  +  (a)  -*  0  as  a  -*  0 
and  condition  (36)  is  satisfied,  for  both  the 
approximations  KM(t)  and  K0(t)  as  may  be  seen  from 
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-Xil X.  a  Kg  y  n  ^  1 
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Fig.  16  —  Steepness  of  the  Transverse  (-)  and  Divergent 
(+)  Waves  in  the  Kelvin  Wake  for  a  Simple  Ship  Bow 
Form,  with  p  =  12°,  y  -  45°  and  F  =  0.5,  in  the 
Michell  Thin-Ship  Approximation  (Left)  and  the  Zeroth- 
Order  Slender-Ship  Approximation  (Right) 


Fig.  17  —  Upper  Bound  for  the  Steepness  of  the 
Divergent  Waves  in  the  Kelvin  Wake  for  a  Simple  Ship 
Bov  Form,  with  p  =  12°  and  y  =  45°,  in  the  Michell 
Thii  ''hip  Approximation  (Left)  and  the  Zeroth-Order 
Slender-Ship  Approximation  (Right) 


equations  (45)  and  (67b).  Figure  16  also  shows  a  sharp 
peak  in  the  steepness  of  the  divergent  waves  at  the  value 
of  a  equal  to  tan/l/(2  +  tan2/)). 

The  steepness  s  +  (a)  of  the  divergent  waves  is  given 
by  equations  (32)  and  (33).  An  upper  bound  for  the 
functions  s  +  (o)  may  be  obtained  by  using  an  upper 
bound  for  the  function  |K(t  +  )|  in  equation  (32).  Equation 
(45)  yields  the  following  upper  bound  for  the  function 
|KM(t)|  corresponding  to  the  Michell  thin-ship 
approximation:  |KM(t)|  <  4/)0/t3  as  t  -*  ».  An  upper 
bound  for  the  function  |K^(t)|  is  given  by  equation  (66), 
where  the  upper  bound  defined  by  equation  (68)  is  used 
for  the  term  |N|.  These  upper  bounds  for  the  functions 


|KM(t)|  and  IKqO)!  can  be  expressed  in  terms  of  a  by 
using  equation  (12).  The  corresponding  upper  bounds  for 
the  steepness  function  s+(o)  are  depicted  in  figure  17. 
Comparison  of  figures  16  and  17  shows  that  the  upper 
bound  for  the  steepness  of  the  divergent  waves  in  the 
Kelvin  wake  depicted  in  figure  17  is  satisfactory  for  all 
values  of  a  for  the  Michell  approximation  KM,  whereas 
that  corresponding  to  the  slender-ship  approximation  Kg  is 
satisfactory  for  values  of  a  smaller  than  approximately 
half  the  value  tan/J/(2  +  tan2/)).  In  both  cases,  the  upper 
bounds  for  the  function  s+(o)  depicted  in  figure  17  are 
satisfactory  for  the  range  of  small  values  of  a  for  which 
the  steepness  s + (a)  is  large.  It  is  noteworthy  that  these 
upper  bounds  for  the  steepness  of  the  divergent  waves  are 
valid  for  all  Froude  numbers,  since  equations  (12),  (32) 
and  (33),  and  the  upper  bounds  for  |KM(t)|  and  |K0(t)|  do 
not  involve  the  Froude  number. 

Figure  18  depicts  the  steepness  function  s+(o)  of 
the  divergent  waves  corresponding  to  the  slender-ship 
approximation  Kfl(t)  for  the  simple  bow  shape  considered 
previously  in  the  two  cases  when  the  maximum  flare  angle 
y  is  taken  equal  to  0  and  45°  (at  the  top  and  bottom 
halves  of  the  figure,  respectively)  and  for  two  values  of  the 
Froude  number,  namely  for  F  =  0.8  and  0.3  (on  the  right 
and  left  halves  of  the  figure,  respectively).  Comparison  of 
the  right  and  left  halves  of  figure  18  shows  no  appreciable 
difference  between  the  values  of  the  steepness  function 
s+(o)  for  the  range  of  small  values  of  a  for  which  the 
steepness  is  large,  in  agreement  with  the  previously-noted 
result  that  the  upper  bound  for  s+(o)  depicted  on  the 
right  side  of  figure  17  is  independent  of  the  Froude 


Wake  for  a  Simple  Ship  Bow  Form  with  p  =  12°,  y  =  0 
(Topi  and  45°  (Bottom),  and  F  -  0.3  (Left)  and  0.8 
(Right)  in  the  Zeroth-Order  Slender-Ship  Approximation 
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Fig.  19  —  Amplitude  of  the  Peak  Value  of  the  Steepness 
of  the  Divergent  Waves  in  the  Kelvin  Wake  ‘or  a  Simple 
Ship  Bow  Form  with  /)  =  12°  and  y  =  45°  in  the 
Zeroth-Order  Slender-Ship  Approximation 

number.  Comparison  of  the  top  and  bottom  halves  of 
figure  18,  corresponding  to  y  =  0  and  45°  as  was  already 
noted,  shows  appreciable  differences,  especially  for  F  = 
0.3  for  which  the  function  s+(o)  exhibits  a  very 
pronounced  peak  in  the  case  y  =  45°.  No  such  peak  is 
apparent  in  the  top  half  of  the  figure  for  y  =  0,  that  is 
in  the  case  when  the  hull  intersects  the  free  surface 
orthogonally. 

The  right  half  of  Figure  16  and  the  bottom  half  of 
figure  18  show  that  the  magnitude  of  the  peak  in  the 
steepness  function  s  +  (o)  for  a  =  tan/?/(2  +  tan2/3)  strongly 
depends  on  the  value  of  the  Froude  number.  Specifically, 
the  peak  is  very  pronounced  in  figure  18  for  F  =  0.3, 
fairly  pronounced  in  figure  16  for  F  -  0.5,  and  almost 
nonapparent  in  figure  18  for  F  =  0.8.  The  magnitude  of 
the  peak  in  the  steepness  function,  that  is  the  value  of  the 
function  s+(o)  for  a  =  tan/?/(2  +  tan2/3),  is  represented  in 
figure  19  as  a  function  of  the  Froude  number,  which  is 
based  on  the  length  of  the  bow  region.  This  figure  shows 


that  the  magnitude  of  the  peak  increases  very  rapidly  as 
the  Froude  number  decreases  below  a  certain  threshold 
value  in  the  vicinity  of  F  =  0.3. 

Figure  20  depicts  the  boundary  of  the  Kelvin  wake, 
which  corresponds  to  a  =  -y/x  =  l/23/2  (that  is,  an 
angle  equal  to  approximately  19°28'),  the  line  a  = 
tan0/(2+tanJ/l)  ~  p/2  (that  is,  an  angle  equal  to 
approximately  6°)  atong  which  the  steepness  of  the 
divergent  waves  has  a  peak,  and  the  lines  along  which  the 
steepness  of  the  divergent  waves  is  equal  to  1/20,  1/15, 
and  1/7  (shown  as  a  chain  line  close  to  the  track  of  the 
ship).  The  latter  three  lines  were  determined  by  using  the 
upper  bound  for  the  steepness  function  s+(a)  that  was 
determined  previously  and  depicted  on  the  right  half  of 
figure  17.  The  four  lines  inside  the  Kelvin  wake  shown  in 
figure  20  correspond  to  the  zeroth-order  slender-ship 
approximation  K0(t)  for  the  simple  bow  shape  considered 
previously,  with  /}  =  12°,  y  =  45°  and  d  =  0.1.  The 
three  lines  along  which  the  divergent  waves  are  steep  lie 
much  closer  to  the  track  of  the  ship  than  the  lines  a  cy  6° 
along  which  the  steepness  of  the  divergent  waves  exhibits 
a  peak. 


Fig.  20  —  The  Kelvin  Cusp  Line,  the  Line  Along  Which 
the  Amplitude  of  the  Divergent  Waves  Exhibits  a  Peak, 
and  the  Three  Lines  Along  Which  the  Steepness  of  the 
Divergent  Waves  is  Equal  to  1/20,  1/15  and  1/7  (Chain 
Line  Close  to  the  Track  of  the  Ship)  for  a  Simple  Ship 
Bow  Form,  with  ft  =  12°  and  y  =  45°,  in  the  Zeroth- 
Order  Slender-Ship  Approximation 

The  three  lines  along  which  the  steepness  of  the 
divergent  waves  is  equal  to  1/20,  1/15  and  1/7,  which  are 
depicted  in  figure  20  for  0  >  x  >  - 100  and  0  <  y  <  40, 
are  represented  again  in  figure  21  at  a  distorted  scale 
where  0  >  x  >  -  300  and  0  <  y  <  7.  The  corresponding 
constant-steepness  lines  predicted  by  the  Michell  thin  ship 
approximation  are  also  shown  in  figure  21  for 
comparison.  The  latter  lines  were  determined  from 
equations  (32)  and  (33)  and  the  upper  bound  |KM(t)|  < 
4/?0/'t3  given  by  equation  (45).  Figure  21  shows  significant 
differences  between,  the  constant  steepness  lines  predicted 
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Fig.  21  —  Lines  Along  Which  the  Steepness  of  the 
Divergent  Waves  is  Equal  to  1/20,  1/15  and  1/7  for  a 
Simple  Ship  Bow  Form,  with  p  =  12°  and  y  =  45°,  in 

the  Michell  Thin-Ship  Approximation  ( - )  and  the 

Zeroth-Order  Slender-Ship  Approximation  ( - ) 


Fig.  22  —  Lines  Along  Which  the  Steepness  of  the 
Divergent  Waves  is  Equal  to  1/15  for  a  Simple  Ship  Bow 
Form  with  P  =  8°,  12°,  16°  and  y  =  45°  (Top)  and 
y  =  0,  25°,  45°  and  p  =  12°  (Bottom)  in  the  Zeroth- 
Order  Slettcler-Ship  Approximation 

by  the  Michell  thin-ship  approximation  KM(t)  and  the 
zeroth-order  slender-ship  approximation  K0(t).  This  figure 
strongly  suggests  the  need  for  performing  additional 
calculations  based  on  a  more  realistic  mathematical  model 
than  the  simple  thin-ship  and  slender-ship  approximations 
used  in  this  study. 

Finally,  the  effect  of  the  entrance  angle  fl  and  of 
the  flare  angle  y  on  the  steepness  of  the  divergent  waves  is 
illustrated  in  figure  22.  Specifically,  this  figure  depicts  the 
lines  along  which  the  steepness  of  the  divergent  waves,  as 
predicted  by  the  slender-ship  approximation  K^t),  is  equal 
to  1/15  for  p  =  8°,  12°,  16°  and  y  =  45°,  in  the  top 
half  of  the  figure,  and  for  y  =  0,  25°,  45°  and  /J  =  12° 
in  the  bottom  half  of  the  figure.  This  figure  shows  that 


the  short  divergent  waves  in  the  Kelvin  wake  become 
steeper  as  the  entrance  angle  p  increases  and  as  the  flare 
angle  y  decreases.  More  generally,  figure  22  shows  that  the 
short  divergent  waves  in  the  Kelvin  wake  are  strongly 
influenced  by  the  hull  shape,  and  it  therefore  suggests  the 
need  for  performing  additional  calculations  in  whioh 
systematic  variations  in  hull  shape  are  considered. 

SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

Asymptotic  expressions  for  determining  the  velocity 
potential  and  its  derivatives  at  a  sufficiently-large  distance 
behind  a  ship  advancing  at  constant  speed  in  calm  water 
are  given  by  equations  (23),  (23a-e),  (24a, b),  (25a-d),  (20a- 
e)  and  (26a-d).  The  far-field  asymptotic  approximation 
(23)  is  uniformly  valid  in  the  vicinity  of  the  track  of  the 
ship  a  =  0  if  condition  (27)  is  satisfied.  For  the  simple 
bow  shape  considered  in  this  study,  condition  (27)  is 
satisfied  for  k  =0  and  1,  corresponding  to  the  potential 
4  and  the  frce-surfacc  elevation  4„,  when  the  far-field 
wave-amplitude  function  K(t)  is  approximated  by  the 
Michell  thin-ship  approximation  KM(t)  or  the  zeroth-order 
slender-ship  approximation  K^t).  However,  condition  (27) 
is  not  satisfied  for  k  >  4  and  k  >  2  for  the 
approximations  KM  and  K0,  respectively. 

The  asymptotic  approximations  used  in  this  study 
provide  simple  explicit  analytical  expressions  for 
determining  the  velocity  potential  and  its  derivatives  for 
large  values  of  -  Xg/U2,  that  is  in  the  far  field,  in  terms 
of  the  far-field  wave-amplitude  function.  However,  for 
small  and  intermediate  values  of  Xg/U2,  these  asymptotic 
approximations  are  not  useful,  and  the  integrals  (18)  must 
be  evaluated  numerically.  For  intermediate  values  of 
Xg/U2,  the  exponential  function  E+(tix)  +  E_(tix)  in  the 
integrands  of  the  integrals  (1),  (3a, b)  and  (5a, b)  oscillates 
fairly  rapidly,  as  may  be  seen  from  figure  1.  Accurate  and 
efficient  integration  rules  suited  to  oscillatory  integrands  of 
the  type  depicted  in  figure  1  must  be  used.  For  small 
values  of  Xg/U2,  on  the  other  hand,  the  oscillations  of 
the  exponential  function  E  +  (tix)  +  E_(t;x)  are  not 
significantly  more  rapid  than  the  oscillations  of  the  fat- 
field  wave-amplitude  function  K(t)  which  also  appears  in 
the  integrands  of  the  wave  integrals  (18),  so  that  a 
different  integration  rule  is  required. 

The  amplitude,  a+,  of  the  divergent  waves  in  the 
Kelvin  wake  vanishes  at  the  track  of  the  ship  if  condition 
(37)  is  satisfied.  However,  it  is  well  known  that  the 
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wavelength,  A+,  of  the  divergent  waves  also  vanishes  at 
the  track  of  the  ship,  as  may  be  seen  from  figures  6  and 
7.  Therefore,  the  divergent  waves  can  theoretically  become 
infinitely  steep  at  the  track  of  the  ship.  More  precisely, 
the  steepness,  s+  =  a + /A  + ,  of  the  divergent  waves  is 
unbounded  at  the  track  of  the  ship  if  condition  (36)  is 
satisfied.  Both  conditions  (37)  and  (36)  can  be  satisfied 
simultaneously  if  condition  (39)  is  satisfied,  where  the  far- 
field  wave-amplitude  function  is  of  order  I/tM  as  t  -*  ». 
Conditions  (37)  and  (36),  and  consequently  also  condition 
(39),  are  satisfied  in  the  cases  of  the  thin-ship  and  the 
slender-ship  approximations  KM  and  K^  for  the  simple 
ship  form  considered  in  the  study,  as  may  be  verified 
from  figures  15  and  16  where  a+(o)  -*  0  and  s+(o)  -*■  » 
as  a  •*  0. 

Infinitcly-steep  water  waves  cannot  exist  in  reality. 
Indeed,  there  exists  a  theoretical  upper  bound  for  the 
steepness  of  water  waves  in  deep  water,  which  is 
approximately  equal  to  1/7.  Condition  (39)  thus  suggests 
that  no  divergent  waves  can  exist  within  a  certain  region 
in  the  vicinity  of  the  track  of  the  ship,  and  that  the  Kelvin 
wake  contains  three  distinct  regions:  (i)  an  inner  region 
adjacent  to  the  track  of  the  ship  where  only  transverse 
waves  can  exist,  (ii)  an  outer  region  where  both  transverse 
and  divergent  waves  are  present,  and  (iii)  an  intermediate 
region  at  the  boundary  between  the  inner  and  outer 
regions  where  steep  short  divergent  waves,  as  well  as 
transverse  waves,  can  be  found. 

Surface-tension  and  nonlinear  effects  have  been 
ignored  in  the  analysis  presented  in  this  study.  This  linear 
no-surface-tension  analysis  predicts  extremely  short  and 
steep  waves  in  the  vicinity  of  the  track  of  the  ship.  Both 
surface-tension  and  nonlinear  effects  therefore  are  liable  to 
be  significant,  and  these  effects  should  be  taken  into 
account.  A  linear  analysis  including  surface-tension  effects 
should  be  performed  first,  since  it  is  evident  from  the 
results  obtained  in  the  present  study  and  from  the  brief 
description  of  the  effects  of  surface  tension  upon  the 
Kelvin  wake  given  in  Lamb  (9,  pp.  468-470]  and 
Wehausen  and  Laitone  [10,  pp.  636-637]  that  the  system 
of  divergent  waves  in  the  vicinity  of  the  track  of  the  ship 
is  likely  to  be  profoundly  affected  by  surface  tension. 

It  was  previously  found  by  Scragg  (51  that,  for  a 
ship  bow  form  with  a  large  flare  angle,  the  zeroth-order 
slender-ship  approximation  K0(t)  predicts  a  sharp  peak  in 
the  value  of  the  amplitude  of  the  divergent  waves  at  a 


value  of  a  equal  to  approximately  half  the  entrance  angle 
p.  This  finding  of  Scragg  has  been  verified  in  this  study, 
as  may  be  seen  from  figures  15,  16  and  18.  Furthermore, 
the  magnitude  of  the  steepness  of  the  divergent  waves  has 
been  found  to  increase  very  rapidly  as  the  Froude  number 
decreases  below  a  certain  threshold  value,  as  is  shown  in 
figure  19. 

The  line  along  which  the  steepness  s+(o)  of  the 
divergent  waves  has  a  peak  and  the  lines  along  which 
s+(e)  takes  the  large  values  1/7,  1/15  and  1/20  have  been 
determined,  for  a  simple  ship  form,  on  the  basis  of  both 
the  zeroth-order  slender-ship  approximation  K^ft)  and  the 
Michell  thin-ship  approximation  KM(t).  Figure  20, 
corresponding  to  the  slender-ship  approximation  1^(0, 
shows  that  these  lines  are  well  inside  the  Kelvin  angle,  and 
that  the  large-steepness  lines  are  much  closer  to  the  track 
of  the  ship  than  the  line  corresponding  to  the  peak  in  the 
steepness  of  the  divergent  waves. 

The  lines,  depicted  in  figures  20  and  21,  along 
which  the  steepness  of  the  divergent  waves  takes  large 
constant  values  are  independent  of  the  value  of  the 
Froude  number,  but  they  strongly  depend  on  the  hull 
shape,  as  may  be  seen  from  figure  22  where  "constant- 
steepness  lines’*  corresponding  to  several  values  of  the 
entrance  angle  p  and  of  the  flare  angle  y  arc  depicted. 

This  figure  shows  that  the  short  divergent  waves  in  the 
Kelvin  wake  become  steeper  as  the  entrance  angle 
increases  and/or  as  the  flare  angle  decreases. 

It  was  found  that  the  lines,  along  which  the 
steepness  of  the  divergent  waves  takes  large  constant 
values,  predicted  by  the  slender-ship  approximation  K0(t) 
and  the  thin-ship  approximation  KM(t)  are  quite  different 
from  one  another,  as  may  be  seen  from  figure  21.  This 
figure  therefore  indicates  the  need  for  performing 
additional  calculations  based  on  a  more  realistic 
mathematical  model  than  the  simple  thin-ship  and  slender- 
ship  approximations  used  in  this  study.  These  two 
approximations  correspond  to  simple  special  cases  of  the 
Neumann-Kelvin  theory,  which  should  then  be  used.  In 
particular,  it  would  be  useful  to  determine  whether  this 
more  realistic  theory  predicts  that  the  steepness  of  the 
divergent  waves  in  the  Kelvin  wake  exhibits  a  peak  (or 
several  peaks),  as  was  found  by  using  the  slender-ship 
approximation  K0(t)  for  a  ship  bow  form  with  large  flare 
angle.  Figure  21  specifically  demonstrates  the  importance 
of  obtaining  accurate  predictions  of  the  far-field  wave- 


amplitude  function  K(t)  for  large  values  of  t.  Indeed,  the 
short  divergent  waves  in  the  vicinity  of  the  track  of  the 
ship,  that  is  for  small  values  of  a,  are  associated  with  Up¬ 
value  of  the  function  K(t)  for  large  values  of  t,  as  may  be 
seen  from  equations  (23),  (24a),  (25a)  and  (12).  Precise 
knowledge  of  the  asymptotic  behavior  of  the  function  K(t) 
as  t  — •  “  therefore  is  critical. 

The  behavior  of  the  function  K(t)  as  t  —  “  has 
been  determined  analytically  in  this  study  for  the  simple 
case  of  the  thin-ship  and  slender-ship  approximations 
KM(t)  and  K0(t)  for  an  idealized  ship  bow  form.  More 
precisely,  the  asymptotic  behavior  of  the  functions  KM(0 
and  K0(t)  for  the  simple  ship  bow  shape  considered  in  this 
study  is  specified  by  equations  (45)  and  (67a,b).  These 
asymptotic  approximations  show  that  we  have  |KM(t)|  = 
0(l/t3)  and  |K0(t)|  =  0(l/t2)  as  t  —  00 .  They  also  show 
that  the  Michell  thin-ship  approximation  KM(t),  which 
corresponds  to  the  thin-ship  limit  of  the  slender-ship 
approximation  Kfl(t),  is  not  uniformly  valid  in  the  limit  t 
—  oo.  It  is  possible  to  analytically  determine  the  behavior 
of  the  far-field  wave-amplitude  function  K(t)  associated 
with  the  Neumann-Kelvin  theory  for  an  arbitrary  ship 
form,  as  is  shown  in  [12].  Such  an  asymptotic 
approximation  for  the  function  K(t)  for  large  values  of  t 
is  useful  because  it  provides  an  explicit  analytical 
relationship  between  the  hull  form  and  the  Froude 
number,  on  one  hand,  and  the  amplitude  of  the  short 
divergent  waves  in  the  vicinity  of  the  track  of  the  ship,  on 
the  other  hand. 

ACKNOWLEDGMENTS 
This  study  was  funded  by  the  Office  of  Naval 
Technology  sponsored  Exploratory  Development  Surface 
Ship  Wake  Detection  Project  at  the  David  W.  Taylor 
Naval  Ship  R&D  Center.  The  authors  wish  to  thank 
Dr.  Arthur  Reed  and  Mr.  Seth  Hawkins  for  their  interest 
in  the  study  and  for  their  useful  comments. 


REFERENCES 

1.  Fu,  Lee-Lueng  and  Benjamin  Holt,  “Seasat 
Views  Oceans  and  Sea  Ice  with  Synthetic  Aperture 
Radar”.  JPL  Publication  81-120,  15  Feb  1982. 

2.  McDonough,  Robert  N.,  Barry  E.  RafFand 
Joyce  L.  Kerr,  "Image  Formation  from  Spaceborne 
Synthetic  Aperture  Radar  Signals”,  Johns  Hopkins  APL 
Technical  Digest,  Vol.  6,  No.  4,  Oct-Dec  1985, 

pp.  300-312. 

3.  Swanson,  Claude  V.,  “Radar  Observability  of 
Ship  Wakes”,  Applied  Physics  Technology  Report  No.  1, 
May  1984,  109  pp. 

4.  Case,  K.M.,  et  al,  “Seasat  Report”,  MITRE 
Corporation,  Report  JSR-83-203,  March  1984. 

5.  Scragg,  Carl  A.,  “A  Numerical  Investigation  of 
the  Kelvin  Wake  Generated  by  a  Destroyer  Hull  Form”, 
Science  Applications,  Report  No.  SA1-83/I216,  Oct  1983, 
pp.  46. 

6.  Cooper,  A.L.,  “Interactions  Between  Ocean 
Surface  Waves  and  Currents”,  Naval  Research 
Laboratory,  NRL  Memorandum  Report  5755,  April  1986, 

pp.  18. 

7.  Noblesse,  F.,  “A  Slender-Ship  Theory  of  Wave 
Resistance”,  Journal  of  Ship  Research,  Vol.  27,  No.  1, 
March  1983,  pp.  13-33. 

8.  Sharma,  S.D.,  “Some  Results  Concerning  the 
Wavemaking  of  a  Thin  Ship”,  Journal  of  Ship  Research, 
Vol.  13,  1969,  pp.  72-81. 

9.  Lamb,  H.,  “Hydrodynanics”,  Dover 
Publications,  New  York,  1879,  pp.  '.'38. 

10.  Wehausen,  J.V.  and  E.V.  Laitonc,  “Surface 
Waves”,  in  Encyclopedia  of  Physics,  Sptingcr-Verlag, 
Berlin,  Vol.  IX,  1960,  pp.  446-778. 

11.  Urscll,  F.,  “On  Kelvin’s  Ship-Wave  Pattern”, 
Journal  of  Fluid  Mechanics,  Vol.  8,  1960,  pp.  418-431. 

12.  Noblesse,  F.,  "Analytical  Approximations  for 
Steady  Ship  Waves  at  Low  Froude  Numbers",  21st  ATTC 
Resistance  and  Flow  Committee,  Washington,  Aug.  1986. 


19 


INITIAL  DISTRIBUTION 


Copies 

Copies 

1 

US  Army  Waterways 

1 

NUSC  NPT 

Experiment  Station 
Research  Center  Lib 

1 

NUSC  NLONLAB 

3 

CHONR/Code  432 

1 

NAVSHIPYD  BREM/Lib 

1 

CONR/Boston 

1 

NAVSHIPYD  CHASN/Lib 

1 

CONR/Chicago 

1 

NAVSHIPYD  MARE/Lib 

1  . 

CONR/Pasadena 

1 

NAVSHIPYD  NORVA/Lib 

1 

NAVPGSCOL 

1 

NAVSHIPYD  PEARL/Lib 

1 

NROTC  &  NAVADMINU 

1 

NAVSHIPYD  PHILA 

1 

NAVWARCOL 

12 

DTIC 

5 

NRL 

1 

AFFDL/FDDS/J.  Olsen 

1  Lib 

4  5841 

2 

AFFDL/FYS 

1  Dale  Cooley 

.i  10 

NAVSEA 

1  S.J.  Pollock 

1  SEA  003 

1  SEA  0P  & 

1  SEA  31B 

1  SEA  55W 

2  SEA  55W3 

2  SEA  55W32 

2 

COGARD 

1  COM  (E),  STA  5-2 

1  Div  of  Merchant  Marine 
Safety 

LC/Sci  &  Tech  Div 

1  SEA  55W33 

1  SEA  56X1 

1 

MARAD/Adv  Ship  Prog  Office 

1 

1 

1 

NAVFACENGCOM 

NAVOCEANO/Lib 

1 

MMA/Tech  Lib 

NASA  AMES  RESEARCH  CENTER/ 

1 

1 

NADC 

R.T.  Medan,  Ms  221-2 

1 

NWC 

4 

NASA  LANGLEY  RESEARCH  CENTER 

1  J.E.  Lemar,  Ms  404 A 

1  Brooks 

1  E.C.  Yates,  Jr.,  Ms  340 

1 

NOSC 

1 

CEL/Code  L31 

1  D.  Bu-‘'  lell 

1 

NSWC/White  Oak/Lib 

1 

NASA/Sci  &  Tech  Info  ’ac  '.ity 

1 

NSWC/Dahlgren/Lib 

1 

NSF/Eng  Div 

Copies 

1  Univ  of  Bridgeport 

Prof.  E.  Uram 
Hech  Eng  Dept 

5  Univ  of  California,  Berkeley 

•  College  of  Eng,  NA  Dept 
1  Lib 

1  J.R.  Paulling 
1  J.V.  Wehausen 
1  W.  Webster 

1  R.  Yeung 

3  CA  Inst  of  Tech 

1  A.J.  Acosta 
1  T.Y.  Wu 
1  Lib 

1  Colorado  State  Univ 

M.  Albertson 
Dept  of  Civ  Eng 

1  Univ  of  Connecticut 

V.  Scottron 
Hyd  Research  Lab 

1  Cornell  Univ 

Grad  School  of  Aero  Eng 

1  Florida  Atlantic  Univ 

Ocean  Eng  Lib 

2  Harvard  Univ /Dept  of  Math 

1  G.  Birkhoff 

1  G.  Carrier 

1  Univ  of  Hawaii /Dr.  Bretschneider 

1  Univ  of  Illinois/Coll  of  Eng 

J.M.  Robertson 
Theoretical  &  Applied  Mech 

4  State  Univ  of  Iowa 

Iowa  Inst  of  Hyd  Research 
1  L.  Land weber 
1  J.  Kennedy 
1  V.ci  Patel 
1  F.  Stern 

1  Kansas  State  Univ 

Eng  Exp  Station/D. A.  Nesmith 


Copies 

1  Lehigh  Univ  Fritz  Lab  Lib 

1  Long  Island  Univ 

Grad  Dept  of  Marine  Sci 
David  Price 

1  Delaware  Univ/Math  Dept 

4  Univ  of  Maryland 

1  Eng  Lib 
1  P 

1  C.L.  S 
1  F.  Buckley 

6  MIT/Dept  of  Ocean  Eng 

1  J.R.  Kerwin 
1  J.N.  Newman 
1  P.  Leehey 
1  M.  Abkowitz 
1  A.T.  Ippen /Hydro  Lab 
1  T.F.  Ogilvie 

5  Univ  of  Michigan/Dept/NAME 

1  Lib 
1  R.  Beck 
1  R.B.  Couch 
1  W.  Vorus 
1  T.E.  Brockett 

5  Univ  of  Minn/St.  Anthony  Falls 

1  R.  Arndt 
1  C.  Farell 
1  J.M.  Ki lien 
1  F.  Schiebe 
1  J.M.  Wetzel 

3  City  College,  Wave  Hill 

1  W.J.  Pierson,  Jr. 

1  A.S.  Peters 
1  J.J.  Stoker 

1  Univ  of  Notre  Dame 

A.F.  Strandhagen 

1  Penn  State  Univ 

Applied  Research  Lab 

2  SAI/Annapolis 

1  N.  Salvesen 
1  C.  von  Kerczek 
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A 


Copies 

3  Southwest  Research  Inst 

1  H.N.  Abramson 
1  G.E.  Transleben,  Jr. 

1  Applied  Mech  Review 

4  Stanford  Univ/Dept  of  Div  Eng 

1  R.L.  Street 

1  B.  Perry 

2  Dept  of  Aero  and  Astro/ 

J.  Ashley 
M.  Van  Dyke 

1  Stanford  Research  Inst /Lib 

2  Stevens  Inst  of  Tech/Davidson  Lab 

1  D.  Savitsky 
1  Lib 

1  Utah  State  Univ/Col  of  Eng 

Roland  W.  Jeppson 

2  Univ  of  Virginia/Aero  Eng  Dept 

1  J.K.  Haviland 
1  Young  Yoo 

2  Webb  Institute 

1  Lib 

1  L.W.  Ward 

1  Worcester  Poly  Inst/Alden 

Research  Lab 

1  Woods  Hole,  Ocean  Eng  Dept 

1  SNAME 

1  Aerojet-General/W.C.  Beckwith 

1  Bethlehem  Steel  Sparrows  Tech  Mgr 

1  Bolt,  Beranek  &  Newman,  MA 

4  Boeing  Company/Aerospace  Group 

1  R.R.  Barberr 
1  W.S.  Rowe 
1  P.E.  Rubbert 
1  G.R.  Saaris 

1  CALSPAN,  Inc.  Applied  Mech  Dept 


Copies 

1  Flow  Research,  Inc. 

1  Eastern  Research  Group 

2  General  Dynamics  Corp- 

1  Convair  Aerospace  Div 
A.M.  Cunningham,  Jr. 

Ms  2851 

1  Electric  Boat  Div 
V.T.  Boatwright,  Jr. 

1  Gibbs  &  Cox,  Inc. 

Tech  Info  Control  Section 

1  Grumman  Aircraft  Eng  Corp 

W.P.  Carl,  Mgr /Grumman  Marine 

2  Tracor  Hydronautics,  Inc. 

1  G.  Miller 
1  A.  Goodman 

4  Lockheed  Aircraft  Corp 

Lockheed  Missiles  &  Space 
1  R.L.  Waid 
1  R.  Lcay 
1  R.  Perkins 
1  R.  Kramer 

1  Marquadt  Corp/F.  Lane 

General  Applied  Sci  Labs 

1  Martin  Marietta  Corp/Rias 

3  McDonnell-Douglas  Corp/Douglas 
Aircraft  Company 
1  Lib 

1  T.  Cebeci 

1  Newport  News  Shipbuilding/Lib 

1  Nielsen,  NA  Rockwell 

1  North  American  Rockwell 
Los  Angeles  Div 
J.R.  Tulinius/Dept  056-015 

2  Northrop  Corp/Aircraft  Div 
1  J.T.  Gallagher 
1  J.R.  Stevens 
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Copies 

Copies 

Code 

Name 

1 

Sperry  Sys  Mgat 

3 

1544 

1 

Robert  Taggart,  Inc. 

.  2 

156 

1 

Tracor 

1 

1561 

2 

VPI 

1 

1562 

1  J.  Schetz 

1  P.  Kaplan 

1 

1563 

1 

1564 

CENTER  DISTRIBUTION 

1 

16 

« 

Copies 

Code  Name 

1 

1603 

1 

012.2 

1 

163 

1 

012.3 

1 

17 

• 

1 

12 

1 

18 

1 

122 

1 

1802 

1 

15 

1 

184 

1 

1502 

4 

1843 

1 

1504 

1 

19 

1 

1505 

1 

194 

1 

1506 

1 

27 

2 

152 

1 

28 

3 

1521 

10 

5211.1 

Reports  Control 

3 

1522 

1 

522.1 

TIC  (C) 

1 

1523 

1 

522.2 

TJC  (A) 

1 

154 

1 

1541 

15 

1542 

3 

1543 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS: 


1.  DTNSRDC  reports,  ft  formal  series,  contain  information  of  permanent  technical 
value.  They  carry  a  consecutive  numerical  identification  regardless  of  their  classification 
or  the  originating  department. 

2.  Departmental  reports,  a  semtformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use. 
They  carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of 
the  originating  department.  Any  distribution  outside  DTNSRDC  must  be  approved  by 
the  head  of  the  originating  department  on  a  case-by-case  basis. 


